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In a recent paper [1[, we introduced "dynamical dark matter," a new framework for dark- matter 
physics, and outlined its underlying theoretical principles and phenomenological possibilities. Un- 
like most traditional approaches to the dark-matter problem which hypothesize the existence of 
one or more stable dark-matter particles, our dynamical dark-matter framework is characterized by 
the fact that the requirement of stability is replaced by a delicate balancing between cosmological 
abundances and lifetimes across a vast ensemble of individual dark-matter components. This setup 
fS) ■ therefore collectively produces a time-varying cosmological dark-matter abundance, and the differ- 

ent dark-matter components can interact and decay throughout the current epoch. While the goal 
c7^ , of our previous paper was to introduce the broad theoretical aspects of this framework, the purpose 

O^ ■ of the current paper is to provide an explicit model of dynamical dark matter and demonstrate that 

pH ■ this model satisfies all collider, astrophysical, and cosmological constraints. The results of this paper 

r—{ , therefore constitute an "existence proof" of the phenomenological viability of our overall dynamical 

dark-matter framework, and demonstrate that dynamical dark matter is indeed a viable alternative 
to the traditional paradigm of dark-matter physics. Dynamical dark matter must therefore be con- 
sidered alongside other approaches to the dark-matter problem, particularly in scenarios involving 
large extra dimensions or string theory in which there exist large numbers of particles which are 
r^ ■ neutral under Standard-Model symmetries. 
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I. INTRODUCTION 

The nature of what constitutes the non-baryonic dark matter in our universe remains one of the most fundamental 
mysteries in particle physics [2\ . The most precise measurements of the relic abundance of this dark matter to date 
are those derived from WMAP data >j] , which yield a value 

ncBuh^ = 0.1131 ±0.0034, (1.1) 

where h ~ 0.72 is the Hubble constant. Beyond this, we know very little about the properties of this dominant 
constituent of the matter density in our universe, save that its interactions with the fields of the Standard Model 
(SM) are extremely weak. One of the reasons why the nature of the dark matter remains so elusive is its apparent 
stability. Observational constraints on the lifetime r^ of any decaying dark-matter candidate x ^^^ quite stringent. 
Indeed, for any particle with a relic abundance ft^^ ^ f^cDM, current limits [4] from cosmic microwave background 
(CMB) measurements, etc., require that 

Tx > 10^^ s . (1.2) 

For this reason, most models of the dark sector posit the existence of a single dark-matter particle (or, in the case 
of certain multi-component dark-matter scenarios (5|, [6|, a small number of such particles) which is either absolutely 
stable (with that stability usually conferred by some additional symmetry, such as R-parity in supersymmetric models, 
KK-parity [3| in universal extra dimensions |8l-[l0|. or T-parity [il[ in little-Higgs theories [ill), or else sufficiently 
long-lived as to satisfy the bound in Eq. (|1.2p . Indeed, the phenomenological consequences of dark-matter decays in 
models with unstable dark-matter candidates [13| can be quite significant. 

Recently, an alternative framework for addressing the dark- matter question has been proposed [l| . In this so-called 
"dynamical dark matter" paradigm, the dark sector comprises not one or merely a few particle species, but rather 
a vast ensemble of different fields (f)i, each of which contributes only a fraction Qi of the total dark-matter relic 
abundance ficDM- None of these fields is presumed to be absolutely stable, and thus a non-zero decay width Ti is 
associated with each field. However, in this framework, the individual relic abundances of the (pi fields are presumed 
to be generated in such a way that the most stable members of that ensemble are the most abundant. By contrast, 
the abundances of the more unstable members are suppressed according to the size of their decay widths. It is this 
balancing between Ti and fli which makes it possible for the phenomenological constraints relating to the effects of 
dark-matter decays to be satisfied. 

In Ref. [l| , we focused on the model- independent aspects of our dynamical dark-matter framework, discussing its 
broad theoretical properties, without any detailed phenomenological analysis or comparison with data. By contrast, 
in this work, as a "proof of concept," we provide an explicit model of dynamical dark matter. In this model, the fields 
which collectively constitute the dynamical dark-matter ensemble are the KK excitations of a light axion-like field 
propagating in the bulk of a spacetime with one or more large, flat extra dimensions. In this model, the fields of the 
SM, as well as the gauge fields associated with some additional, non-Abelian gauge group G which confines at a scale 
Ag, are taken to be localized on a four-dimensional subspace of that bulk. The axion field is assumed to couple to the 
gauge fields of G (and also potentially to one or more of the SM fields) via non-renormalizable operators suppressed by 
some effective, four-dimensional cutoff scale fx- We shall demonstrate that within this setup, the resulting ensemble 
of axion KK modes naturally satisfies all applicable observational constraints on dark-matter decays — even if the 
stability of this "dark tower" is entirely unprotected. 

Another advantage of this particular dynamical dark-matter model is that it is not only phenomenologically viable, 
but also theoretically well motivated. In any theory in which the SM fields reside on a brane, the KK excitations of 
any bulk field are, from the point of view of the four-dimensional theory, massive particles neutral under the SM gauge 
group. Thus, were it not for the lack of a stabilizing symmetry, any of these particles would be a natural candidate 
for dark matter. However, our results demonstrate that a lack of stability is not an insurmountable impediment to 
such fields serving as dark matter collectively, rather than individually. This is an exciting prospect, for it provides 
a novel way of addressing the dark-matter question in theories with extra dimensions. Furthermore, our model also 
demonstrates that realizing a viable dynamical dark-matter ensemble does not require an overly complicated dark 
sector, a large number of independent mass scales, or an excessive degree of fine-tuning. Indeed, the model presented 



here involves only three independent physical scales: the effective four-dimensional cutoff scale fx, the confinement 
scale A(3, and the compactification scale M^,. Together, these three scales determine the mass spectrum and decay 
properties of the entire ensemble. 

The outline of this paper is as follows. In Sect.|lll we briefly review the formalism for discussing axions and axion-like 
fields, beginning with the standard, four-dimensional case and then moving on to the generalized, five-dimensional 
bulk-axion case. In Sect. IIIH we calculate the decay widths of the KK modes of such a bulk axion-like field and 
investigate how these decay widths scale with the mass of the mode. In the process, we show that the decays of the 
lighter modes to SM fields experience a natural suppression, but that such decays nevertheless dominate over decays 
to other, lighter bulk fields in the theory. In Sect. lIVi we examine the various mechanisms through which a population 
of axion modes may be generated in the early universe, and demonstrate that the abundances generated for those 
modes by misalignment production are indeed balanced against their decay widths in precisely the manner required 
for dynamical dark matter. In Sect. |Vl we examine the collective properties of the ensemble of axion KK modes. 
We show that such an ensemble can collectively reproduce the observed dark-matter relic density given in Eq. (jl.ip . 
and that it possesses the appropriate equation of state to be regarded as dark matter. In Sect. IVH we summarize 
the experimental, astrophysical, and cosmological constraints on scenarios involving bulk axions in large, flat extra 
dimensions. We demonstrate that these constraints can be satisfied in a model which also simultaneously yields the 
correct total relic abundance — in other words, that our model truly constitutes a viable model of dynamical dark 
matter. Finally, in the Conclusions, we summarize the results of the previous sections and discuss several further 
directions for future investigation. 

As we have indicated, this paper is the second part of a two-part series that began with Ref. [1|. Consequently, 
we shall assume that the reader is familiar with the ideas, notation, and conventions established in Ref. [l[ in what 
follows. 

II. BULK AXIONS AS DYNAMICAL DARK MATTER 

As discussed in the Introduction, the model for dynamical dark matter that we shall consider in this paper is a 
model in which the KK excitations of a bulk axion constitute the dark-matter ensemble. In this section, therefore, 
we briefly review the formalism relevant for describing the dynamics of axions in four or more dimensions. We begin 
with a brief summary of the relevant properties of the four-dimensional QCD axion (more detailed reviews of which 
may be found, e.g., in Refs. [15l4l8j ). and then discuss how this formalism can be generalized to a broader class of 
axions and axion-like fields. Finally, we summarize the formalism for embedding such fields in the bulk in theories 
with extra dimensions. 

A. Axions in Four Dimensions 

The QCD axion emerges as a consequence of the Peccei-Quinn (PQ) mechanism |14| . a mechanism which provides 
an elegant, dynamical solution to the strong CP problem. The strong CP problem arises due to the non-trivial vacuum 
structure of QCD. Specifically, the QCD Lagrangian can in principle contain an additional term 
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where gs is the 517(3) coupling, ^ is a numerical factor of 0(1), G°^ is the field-strength tensor for the gluon field, 
and G" = \e^^paG'''"^ is its dual. The parameter 8 is given by = 0-1- ArgdetM, where G is the strong-interaction 
theta-angle and M is the Cabbibo-Kobayashi-Maskawa (CKM) matrix. In principle, 9 can take any value. However, 
experimental bounds on the electric dipole moment dn of the neutron serve to constrain 8. The most stringent limit 
is currently |d„| < 2.9 x 10~^^ e cm [19|, which translates into a bound 

8 < 0.7 X 10"" . (2.2) 

While there is, in principle, no problem with 8 taking so small a value, there is no particular reason why it should 
be so small. This fine-tuning issue is what is commonly referred to as the strong CP problem. 

In the Peccei-Quinn solution to the strong CP problem, the effective 8-parameter associated with the gluon field 
relaxes to zero dynamically as a consequence of the spontaneous breaking of an anomalous, global U{\) symmetry, 
usually dubbed [/(1)pq, at some high scale /pq. The spontaneous breaking of this C/(1)pq symmetry implies the 



presence of a pseudo-Nambu-Goldstone boson: a real pseudoscalar commonly known as the QCD axion [20|, which 
necessarily interacts with the gluon field via the Lagrangian 

£ 3 ^9,,a5^a+ gJ^'^^^ aG^-G°, , (2.3) 

where a denotes the axion field. Furthermore, this pseudoscalar may also have interactions with the other fields of the 
SM. The presence of the anomalous C/(1)pq symmetry in the high-scale theory determines the effective Lagrangian 
for these interactions (at leading order) to be 

where ipi are the SM fermions and Cy and Ci are dimensionless coefficients. These coefficients depend on the charge 
assignments of the SM fields (and potentially of additional fields in the theory as well) under U{1)pq, and are therefore 
substantially more model-dependent than ^. 

At high temperatures, the axion field is effectively massless, as befits a Nambu-Goldstone boson; however, it acquires 
a small, temperature-dependent mass ma{T) at lower scales due to QCD instanton effects. A number of computations 
of this mass have been performed, and while the results depend to some extent on the assumptions and calculational 
techniques involved, ma{T) is often assumed to have the rough form [2lll23| 

, X 9.^ A?.cD 6(^^^^) forr>AQCD 

ma{T) « -^ ^ x{ \ T ) ^ QCD ^2.5) 

Ay/2-K /PQ f T^ < A 

L 1 for T < Aqcd , 

where Aqcd ~ 250 MeV is the QCD confinement scale, 6 is a numerical coefficient of 0(10^^), and i^ is an 0(1) 
numerical factor. 

The fact that the axion necessarily couples to the gluon field implies that it will also have effective couplings to 
hadrons. The most important such couplings, phenomenologically speaking, are those of the axion to pions and 
nucleons. These couplings take the form [23| 

/:had = -:%^(a^a)[(a^7r+)^-7r" + (5^7r-)^+7r°-2(a^^°)7r+7r-] 

/tt/PQ 

+ ^ {d^a)n^^^fn + ^ {d^a)p^^^^''p + ^%^ {d^a) [k+ [p^n) - ^- [n^p)] . (2.6) 

JPQ JPQ hjPQ 

The precise values for the effective nucleon-nucleon-axion couplings Cap and Can, the nucleon-pion-axion coupling 
Ca-nN, and the axion-pion-pion coupling Ca-n depend on the C/(1)pq charges of the quark fields. For the case of a 
so-called hadronic axion 24], which does not couple directly to the SM quarks, the coefficients for the axion-nucleon- 
nucleon interaction arc 



...-^-^+0.02, C,„ = 0.24 -^+0.15 1-^ 

(1 + z) (1 + z) (1 + z) (1 + z) 



Cap = 0.24 ^_, ^ _^ +0.f5 ^_, , _^ +0.02 , Can = 0.24 ^_, ^ _^ +0.15 ^^ , _^ -F 0.02 , (2.7) 



where z = mu/md ~ 0.56 denotes the ratio of the up-quark to down-quark masses. Similarly, the coefficients for the 
interactions involving pions are given by 

where mTr ~ 135.0 MeV is the neutral pion mass, and f-^ ~ 93 MeV is the pion decay constant. These hadronic 
couplings play an important role in constraining the parameter space of axion models. 

The QCD axion is the prototypical example of a light pseudoscalar field whose mass arises solely due to non- 
perturbative effects associated with instanton dynamics, and whose interactions with the SM fields are highly sup- 
pressed. It is by no means the only example, however. Indeed, a wide variety of additional particles possessing 
these same properties have appeared in the literature in a number of beyond-the-Standard-Model (BSM) contexts, 
and are often generically referred to as axion-like particles (ALPs). One particularly well-motivated example is the 
model-independent axion |25| in string theory. Since axions of this more general sort are, by and large, no less viable 
as dark-matter candidates than the QCD axion, it behooves us to extend our focus to encompass such fields as well. 



For the remainder of this work, then, we will use the term "axion" to refer to any pseudoscalar field whose mass 
is generated by the instanton dynamics associated with an arbitrary non-Abelian gauge group G. This gauge group 
could be the SM SU{?)) color group, as it is for the QCD axion, but alternatively it could be some additional group 
which either resides in a hidden sector, or else confines at a very high scale. As with the QCD axion, any axion 
we consider will be assumed to be a pseudo-Nambu-Goldstone boson associated with the breaking of some global 
symmetry U{l)x at a scale fx by the vacuum expectation value (VEV) of some scalar field. The axion field is 
assumed to couple to the field strength Q^^^ associated with G via a term of the form 

Ant 3 .^aGlJ^^^'^ , (2.9) 

o2ttJx 

where a again denotes the axion field, go is the coupling constant associated with G, Qf^"'^ is the dual of 5^'^'', and ^ 
is a model-dependent coefficient which parameterizes the strength of the effective interaction between the axion and 
the gauge fields. We will also assume that G goes through a confining phase transition at some scale Aq, and that 
a potential analogous to that appearing in Eq. (J2.18I) is thereby generated for a. In other words, this general axion 
couples to G in a manner completely analogous to that in which the QCD axion couples to the SM SU{3). It therefore 
follows that all of the QCD-axion formalism outlined above continues to hold for axions in the broader sense of the 
word, provided one makes the substitutions Aqcd — >■ Aq, g^ — ?► go, /pq — ?► fx, etc., where appropriate. 

There is, however, one crucial physical distinction between axions in general and the QCD axion in specific: for 
general axions, the confinement scale Aq is essentially a free parameter. The properties of such an axion are therefore 
far less constrained than those of a QCD axion, simply because the axion mass is not uniquely determined by fx 
alone. Moreover, the vast majority of the experimental bounds on axions depend crucially on the charge assignments 
of the SM gauge fields under the global U{l)x symmetry. For a generic axion, these charges need not have any 
relationship to the C/(1)pq charge assignments for these fields. An important consequence of this is that a generalized 
axion need not couple directly to the gluon field at leading order. Moreover, other scenarios could be realized in such 
a framework that cannot arise for a QCD axion. For example, one can imagine a purely "photonic" axion which 
couples to the photon field at leading order, but not to the gluon field or to any of the SM fermions. In what follows, 
we will focus on several different concrete coupling scenarios. One of these will be such a photonic axion; another 
will be a "hadronic" axion which couples to the gluon and photon fields, but not directly to any of the SM fermions. 
However, we note that numerous other possibilities exist, and that the laboratory, astrophysical, and cosmological 
constraints on any given model depend sensitively on the couplings between the axion and the fields of the SM. 

It is also worth noting that certain details of any scenario of this sort will depend on the details of the instanton 
dynamics associated with the particular gauge group G in question. The scaling behavior of mx (T) as a function of 
T, for example, may not be identical to the scaling behavior quoted in Eq. (|2.5p for QCD instantons. However, none 
of these details plays a crucial role in the dark-matter phenomenology of the our model. We will therefore assume 
for the remainder of this work that, except for the values of Aq, go, etc., the standard axion results derived in the 
context of QCD-instanton dynamics apply to G-instanton dynamics as well. 

B. Axions in Extra Dimensions 

Having summarized the formalism applicable to a four-dimensional axion, we now consider how the situation 
changes when the axion in question is allowed to propagate in the extra-dimensional bulk of a theory with more than 
four dimensions. As was originally pointed out in Ref. [26[, the dynamics of such an axion is far richer than that 
of a purely four-dimensional axion, due both to the presence of an entire KK tower of axion excitations and to a 
non-trivial mixing between these excitations due to the presence of brane mass terms, which explicitly violate KK 
mode-number conservation. Indeed, as we shall see, it is those KK excitations which will constitute the dark-matter 
ensemble in our model, and it is their mixing which gives this ensemble the appropriate properties to be a viable 
dynamical dark-matter candidate. Of course scenarios involving large extra dimensions have many other attractive 
features as well: they provide a geometric interpretation of the hierarchy between the weak scale and the Planck 
scale 27 J29|, between the weak scale and the grand- unification scale Q, and between the weak scale and the string 
scale 9|, |30|. Moreover, a higher-dimensional axion field can be accommodated quite naturally in such a brane/bulk 
framework. Indeed, while only gravity is required to propagate in the bulk, the propagation of SM-gauge-singlet fields 
there, including axions of all varieties, is, in a sense, almost expected. 

In what follows, we present the setup for a generic axion field in the bulk. This parallels the setup for a QCD axion 
put forth in Ref. |26| . For concreteness, we choose to focus on the case in which the axion is allowed to propagate in a 
single, large extra dimension compactified on a Si/'2.2 orbifold of radius R, while the fields of the SM and the gauge 
fields associated with the additional symmetry group G are confined to the brane located at x^ = 0. However, we 



emphasize that the setup described here can easily be extended to scenarios in which the axion in question is ahowed 
to propagate in multiple extra dimensions, or in which the background geometry is more complicated 31]. 

At scales below the weak scale but above the confinement scale Ag, the effective action for a bulk axion in five 
dimensions takes the form 



'S'off — 



d^x 



2TrR 



da; 5 



-dmad^^a + S{X5) (£brane + -Cint) 



(2.10) 



Here, we have divided the brane- localized terms in the Lagrangian into two parts. The first, £branc, contains the 
terms involving the brane fields alone — both the fields of the SM and any additional fields, including the gauge fields 
associated with the gauge group G. The second, Cmt^ contains the interaction terms involving the brane- localized 
fields and the five-dimensional axion. This second piece is given by 



.^int — 



9U 
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3/2 
X 



aQf^.Q"'" 



E 



„3/2 
J X 



(a^a)V',7'^7'V'. 



glcg 



32n^f 



3/2 
X 



iG^.M"'"' 



32n^f 



3/2 
X 



aF„,„Ft''' 



(2.11) 



where e and gg are the respective couplings for U{1)-em and SU{3) color, and fx is the fundamental five-dimensional 
scale associated with the breaking of U{l)x (the analogue of the Peccei-Quinn scale /pq in Ref. |2q|). 

The first term in iZint is the requisite coupling between the five-dimensional axion a and the gauge fields of G. The 
second term represents the derivative couplings between the five-dimensional axion a and the SM fermion fields ipi, 
with mo del- dependent coefficients Ci that depend on the U{l)x charges of the tpi. The remaining two terms represent 
the interactions between the axion and the gluon and photon fields, the field-strength tensors for which are here 
respectively denoted C^"^" and F'^", with (once again model-dependent) coefficients c^ and Cg. 

The five-dimensional axion field can be represented as a tower of KK excitations via the decomposition 



a{x^, x^) 



V2^^n 

^ Tl— 



^ r„a„(a;^) cos (—^) 



(2.12) 



where the factor 



1 for n = 
\/2 otherwise 



(2.13) 



ensures that the kinetic term for each mode is canonically normalized. Substituting this expression into Eq. (|2.1ip 
and integrating over 2:5, we obtain 



'off 



d^x 



E 



df.and'^an 



32^2/ 



x 



-TnanGlJ^^" + y^ ^r„(9^a„)V',7^7V^ 






IX 



9^8 Cg 
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-r„a„^^.F^" - V{a) 



(2.14) 



where 






(2.15) 



and where the quantity /x, defined by the relation 



/: 



2 _ 

X = 



27rRl 



X : 



(2.16) 



represents the effective four-dimensional J7(l)x-breaking scale. Note that each mode in the KK tower couples to the 
SM fields with a strength inversely proportional to fx- Note also that at these scales, the axion mass-squared matrix 



M 



2 

mn 



d^Vja) 
damdan 



(2.17) 



is purely diagonal. 



The effective action in Eq. (|2.14l) is valid at high scales where T ^ Aq. Around T ~ Aq, however, instanton 
effects give rise to an additional contribution to the effective axion potential. In the low-temperature regime, the full 
potential takes the form 



°° 1 Tl2 



n=0 



2i?2 



9g a 4 

32^2 ^« 



1 - cos ^ ^ r„a„ + Oc 
V-^-^ n=0 



(2.18) 



where Oq is the analogue of the QCD theta-parameter Q. Minimizing the potential yields the vacuum configuration 
(ao) — /x(— ©G + 'r^)/C for ^ ^ 2Z, with (a„) = for all n > 0. This additional potential term modifies the axion 
mass-squared matrix at scales T < Aq to 



M'^ 



Mjn 5. 



9ie^h. „ „„./ i 
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'm'n COb 



X 



^/jf j,^o 



^Tfea/c + Qg 



(2.19) 



where Mc = 1/R is the compactification scale. We see here that the terms originating from the instanton-induced 
potential in Eq. (|2.18p include off-diagonal contributions, which result in mixings among the KK eigenstates. In the 
vicinity of the minimum of T^(a„), the axion mass-squared matrix above takes the form [261] 



7W^ 



i-x 



/ 1 %/2 V2 V2 

V2 2 + y^ 2 2 

^/2 2 2 + V 2 

^/2 2 2 2 + V 



V : 



(2.20) 



where 



y 



mx 



and 



.2 _ 9he^% 



The eigenvalues A^ of this mass-squared matrix are the set of solutions to the transcendental equation 



(2.21) 



^^^^cotf^ = A^ 



(2.22) 



The normalized mass eigenstate ax corresponding to each A may be written as a sum of the KK eigenstates a„: 



a\ 



y^ Uxnttn = y^ 



n=0 



Q y A^ — Ti^y^ 



A-XQn , 



(2.23) 



where A — X/rnx, and where 



A..£ 



1 + A' + TrVy' 



-1/2 



The quantity A^ can be shown to obey the sum rules [26[ 

^A^ = l, E^'^A = 1 



(2.24) 



(2.25) 



which follow directly from the unitarity of Uxn- Upon rewriting Eq. (j2.14p in this mass eigenbasis, we obtain the 
axion effective action at temperatures T < Ac, which, up to 0{a^^/f^), is given by 

5eff = fd^x\Y^(ld,axd'^ax~l\'rnW^+"^^^axF,.F'^'^ + 3-^^^ 
J . \^ ^ 52n^fx 327r2/x 



Y: "-^{d,ax)i^n^j'A) + 5^ E XlX',XrxlAx^Ax^Ax,Ax,ax^ax,ax,ax, 
,. fx J 7687r2/i V..A..A..A. 



(2.26) 
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The quartic axion self-interaction terms shown above originate from the instanton- induced potential in Eq. (|2.18l) . 
Other, higher-order terms not shown may be safely neglected when T <C /x- 

If a non-trivial coupling exists between the bulk axion and the gluon field, effective interactions will also arise 
between the ax and the hadron fields at temperatures below Aqcd- (In the case of the QCD axion, of course, such 
couplings are mandatory.) The Lagrangian which describes these interactions is just the five-dimensional analogue of 
Eq. 



Ca 



>Chad = X^Ax^:^id^ax) (9^^+)^-7r" -t- (a'^7r-)7r+7r" - 2(5'^7r")7r+^- -I- X^Ax^^{d^ax)n-f''j^n 



f-irfx 



Ca 



f. 



X 



+ \^Ax^{d^ax)ln^-i''p + i\^Ax-r^{d^ax) n+jri^^n-^-urP , (2.27) 

JX JttJX '- -' 



where the coefficients CaTr, Can, etc., depend on the details of the theory, and may differ from those given in Eqs. (|2.7p 
and dSH). 

As discussed in Ref. []|, the mass spectrum of the model reduces to a reasonably simple form in certain limiting 
cases which depend on the value of the ratio y defined in Eq. (I2.21[) . The first of these is weakly-mixed regime, in 
which y 3> 1. In this regime, the a\ are all very nearly equivalent to the KK modes a„, with masses A ~ nMc, where 
n is an integer. The extent to which any given value of A differs from Mc is set by the size of the off-diagonal terms in 
Eq. (|2.20p . and in particular, the lightest mass eigenstate axg has a mass Aq « Mc/y = mx- In short, the KK tower 
essentially comprises a single light mode plus a tower of massive KK excitations of that mode. In the extreme limit, 
in which Mc — ^ oo, the theory reduces to an effectively four-dimensional theory with a single light axion whose mass 
is precisely equal to mx, as expected. 

In the opposite limit, in which y <ti I, the situation is markedly different 26]. The lighter mass eigenstates in the 
theory have masses X ^ {n+ l/2)Mc ^ mx (where n = {0, 1, 2, . . .} is an integer), while the heavier mass eigenstates 
have masses A « nM^- The transition region between the two regimes occurs at around 

A ~ Atrans = — TT — , (2.28) 

which corresponds to a value n = tt/j/^. In this regime, the states with masses below this threshold are highly 
mixed, owing to the large, off-diagonal terms in Eq. (|2.20l) proportional to mx- We dub this the "strongly-mixed" 
regime. In this latter regime, as we shall soon see, the full KK tower plays a far larger role in the dark-matter 
phenomenology of a given model than in weakly-mixed scenarios, in which the dark-matter phenomenology is, more 
or less, the phenomenology of the zero mode. Indeed, the bulk-axion scenarios which give rise to dynamical dark- 
matter ensembles in which the full tower contributes significantly tend to be those in which y is small. 

It is evident from Eqs. (I2.26|) and (|2.27p that the parameter combination )? A\l jx plays a critical role in bulk- 
axion dynamics. Indeed, it is this combination which determines the strength of the interaction between a given mass 
eigenstate a\ and any of the SM fields. It turns out to be phenomenologically quite significant that not all a\ couple 
to the fields of the SM with the same strength. We discuss the impact of such a coupling structure on the decay 
properties of the tower states in Sect. IIIIl and summarize its impact on phenomenological constraints in Sect. IVII 
(a more detailed analysis of which can be found in Ref. |35|]). In addition, a plot of how A\ and }?A\ depend on 
A is provided in Ref. [l|. It is worth remarking that the non- universality of the a\ couplings is yet another direct 
consequence of the non-trivial mixing between axion KK modes implied by Eq. (|2.20p . This effect does not arise for 
bulk fields in the absence of such mixing: the couplings of the KK excitations of the graviton to the SM fields in 
theories of this sort, for example, are identical for all modes. 

In summary, our model for dynamical dark matter consists of an axion propagating in the bulk of an extra dimension 
of radius i?, with the SM living on a brane. From the perspective of an observer on the brane, our dynamical dark- 
matter ensemble consists of the KK modes of this bulk axion field. As we have discussed above, our model involves 
three important dimensionful parameters: jx, Mc, and Kq. We have also shown that the physics of this model 
depends crucially on one particular dimensionless combination of these parameters, namely y, which governs the 
extent to which the individual KK modes mix when forming the constituents of our dynamical dark-matter ensemble. 

III. CHARACTERIZING THE CONSTITUENTS: DECAY WIDTHS 

Now that we have reviewed the setup underlying our model for dynamical dark matter, we may begin to assess its 
phenomenological ramifications. As discussed in Ref. 1], the essence of the dynamical dark-matter framework lies in 
the balance between the decay widths and relic abundances of the fields which contribute to ficDM- Therefore, our 



principal aim must be to evaluate the decay widths Tx and relic abundances il^ of the fields ax which our dark-matter 
ensemble comprises, and examine how these two quantities scale with A. In this section, we focus on decays: we 
calculate the partial-width contribution associated with each of the potentially relevant decay channels for a generic 
ax and assess how the total width Tx depends on the dimensionful parameters of the model, namely fx, M^, and Aq. 
In the subsequent section, we focus on abundances. 

A. Decays to Standard-Model States 

We begin our discussion of axion decays by computing the decay widths of the ax directly to SM states on the 
brane. The first step is to derive Feynman rules for the relevant interactions, which can be obtained directly from the 
interaction terms given in Eq. (j2.26|) and (|2.27|) . For those ax with A below a few GeV, the relevant vertices are: 
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-X^Axka ■ (fc_ + fc+ - 2fco) 



Here, and throughout the rest of this work, the symbol Mp represents the effective, four-dimensional reduced Planck 
mass. Since the coupling of each ax to the SM fields is suppressed by fx, the contribution to the total decay width 
Tx of each ax which comes from decays to SM fields will be suppressed by fx'^. At low temperatures, the most 
relevant decay processes (depending, of course, on the precise values of Ci and c-y in any given model) will be ax — >■ 77, 
ax -^ e+e~, and ax — >■ i^ii^i, where i — {1, 2, 3} labels the three light neutrino mass eigenstates. 

The one decay channel which is kinematically accessible for all axion modes, regardless of their mass (provided c^ 
is non- vanishing) , is a — >■ 77. The partial decay width of an axion mass eigenstate ax into a pair of photons is 



T{ax -^ 77) 



c^g^A^ 
2567r3/| 



{X'Axf - G-,{\'Axy 



A^ 



(3.1) 



X 



where a = e'^/An, and where we have defined the quantity G^ = c'^a'^ /256n^. Note that r{ax -^ 77) includes 

an overall factor (X^Ax)^- This factor is a direct consequence of mixing in the axion mass matrix, and it appears 
universally in all partial-width expressions for axion decays to states on the brane. For those modes with masses 
A > TTnix/Mc, this factor is 0(1); however, for those modes with A < nm^/Mc, this factor can be much smaller. We 
therefore see that the mixing factor (X^Ax)'^ suppresses the decay widths of the lighter ax, while leaving the widths 
of the heavier ax unsuppressed. This decay-width suppression for the light modes plays a crucial role in bulk-brane 
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models of dynamical dark matter, as discussed in Ref. |l|: since these light modes also turn out to have the largest 
relic abundances, their decays are generally the most dangerous from a phenomenological perspective. 

As A increases, additional decay channels open up in which a given a^ decays to a light fermion-antifermion pair, 
provided that direct couplings exist between the axion and the fermionic species ipi in question. The partial width 
for any decay of this sort is given by 



T{ax -^ V»V'») 






A2 



1/2 



(3.2) 



where m^ is the mass of the fermion in question. Note that if the five-dimensional axion field couples to the SM 
neutrinos, the Tx will have a non-trivial dependence on the neutrino mass spectrum. While the precise masses of the 
three neutrino species are as yet unknown, measurements of the solar and atmospheric squared-mass splittings place 
lower limits on two of the three m^. . The current best-fit values for these splittings are |33 | 



lAmil 



7-59J:J]:2? X 10"^ eV^ 
' eV^ 
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(3.3) 



In what follows we shall assume a normal hierarchy, and we will take the mass of the lightest neutrino to be vanishingly 



small. In this case, the masses of the heavier two neutrinos are m,. 



8.7 X 10"'^ eV and m. 



are comparable to the lower bound [33|| 



Mr > 3.9 X 10"'^ eV 



4.9 X 10"^ eV, which 



(3.4) 



on the compactification scale Mc from modified-gravity experiments. Therefore, if Ale lies only slightly above this 
bound, the masses A of certain light ax will be comparable to m^^ 3 ■ The partial width for the decays of those ax 
to neutrinos can therefore in principle be quite large compared to their partial widths for decays to photons, as can 
be seen by comparing Eqs. (|3.ip and (|3.2p . Note that if c,^. = for all i, the decay width loses all sensitivity to the 
neutrino mass spectrum. This is indeed the case for the photonic and hadronic axions which will serve as our primary 
examples in what follows. 

As A increases still further, decays of the ax to hadrons become kinematically accessible — provided, of course, 
that either Cg ^0, or else that Cq^ ^ for some quark species qt. The lowest such threshold is that for decays of the 
form a — >■ 7r"'"7r~7r'^, which are kinematically allowed whenever A > 277i.n.± -I- m'^. The relevant interaction vertex is 
that appearing in the top line of Eq. (j2.27p , and the corresponding contribution to the decay width of ax from this 
three-body decay takes the form 



rc„ , _+_- J^ C!^^{\ Ax) , , 

1 (ax -> TT^TT TT ) =: , , , ,0 ?■> ^y-^) 



where 1(A) denotes the phase-space integral 
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(3.6) 



with integration variable m^2 = (fc+ + fc-)^. For A ^ 2m^± + m^o, this expression takes the asymptotic form 



T{ax -)■ TT+TT-Tr") w (2.07 x 10"^) q2 
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(3.7) 



In practice, this asymptotic expression is a good approximation for T{ax — >■ tt+tt tt^) as long as A roughly exceeds a 
few GeV. 

For even larger values of A, decays to nucleons become kinematically accessible. In the present treatment, however, 
any ax with masses this large will not play a significant role in the phenomenology of such "dark-tower" scenarios, 
nor will they have a significant impact on the observational and experimental constraints on such scenarios. Indeed, 
as was observed in Ref. JSJ] for the related case of KK-graviton decays, such modes are innocuous precisely because of 
the large contributions to their decay widths from hadronic decays. The axion case under consideration here differs 
qualitatively from this KK-graviton case only in that the decay width F^ of each ax contains an additional factor 
(X^Ax)^- However, as we shall see in Sect. I VII the quantity iTm^/Mc is never much larger than a few GeV in realistic 
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axion models of dynamical dark matter. Consequently, this factor will be 0(1) for modes with A roughly exceeding a 
few GeV — those modes for which decays to nucleons are kinematically allowed — and the partial widths for those 
decays will therefore be unsuppressed. We therefore refrain from explicitly calculating the partial-width contribution 
from such decays, and simply acknowledge that any ax with A roughly exceeding a few GeV will decay quite early — 
i.e., before the big-bang-nucleosynthesis (BBN) epoch. 

Thus far, we have obtained partial-width expressions for all of the relevant decay channels for those a\ with A less 
than roughly a few GeV directly into final states involving SM fields alone. The total width for the decays of a\ into 
this set of final states is obtained by combining the expressions in Eqs. p.ip . (I3.2p . and (|3.5p : 
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where the Heaviside functions 0(A — 2m^) and 0(A — 2m^± — m^o) enforce that only kinematically accessible decays 
contribute in the sum. 

B. Intra-Ensemble Decays 

Up to this point, we have only been considering the decay- width contributions from decays of a given ax directly 
to the fields of the SM. We have yet to address the issue of intra- ensemble decays — i.e., decays in which a given 
state in the dynamical dark-matter ensemble decays to a final state containing one or more other, lighter states in 
that ensemble. Indeed, as discussed in Ref. 'V\, such decays can have a significant impact on the phenomenology of a 
given dynamical dark-matter model: not only will they alter the individual relic abundances fli of the particles (j)i in a 
given ensemble, but they will also alter the phase-space distributions fi{pi, t) of those particles, potentially generating 
a sizable population of (pi with relativistic three-momenta pi. For these reasons, it is crucial to assess whether such 
decays occur at a substantial rate in the bulk-axion model presented here, or whether the net rate for these decays is 
negligible, in which case the quantity T\ given in Eq. (|3.8p truly represents the total decay width of a given a\ with 
A less than a few GeV. 

In the model under discussion here, the dark sector properly comprises KK axions, KK gravitons, and graviscalars. A 
complete description of the dynamics of the ensemble would therefore involve solving the coupled system of Boltzmann 
equations which describes the evolution of the respective phase-space distributions fxipx, t), fniPn, t), and fsiPs, t) for 
the various axion modes ax, KK gravitons G/Z , and graviscalars ips, as discussed in the Appendix of Ref. jli]. In the 
present work, our aim will not be to solve the Boltzmann equations in complete generality, but rather to demonstrate 
that the effects of intra-ensemble decays on the abundances and phase-space distributions of the (j)i are sufficiently 
small that they may be safely neglected for any otherwise phenomenologically reasonable choice of model parameters. 
Some of these effects — for example, the depletion of fix for a given a^ due to intra-ensemble decays — depend only 
on the net contribution T\ to the width of a given a^, obtained by summing over the partial widths from all such 
decays. On the other hand, certain other effects, such as the increase in the abundances f^A and alteration of the 
phase-space distributions fx {px , t) of the lighter ax due to the decays of the heavier ax , depend on these partial widths 
in a different manner. Determining the precise magnitude of these effects therefore requires a more thorough analysis 
of the Boltzmann equations. Nevertheless, T-^ can still be useful as a rough rubric for assessing whether or not the 
effects in question are likely to be significant. In this work, then, we simply demonstrate that T^ for any given 
ax is sufficiently small in comparison with the result in Eq. p.8p in otherwise phenomenologically reasonable regions 
of model parameter space that its effect on total decay widths may safely be neglected. This result provides a good 
initial indication that the additional effects mentioned above are also unimportant. A more rigorous justification for 
neglecting these effects will appear in Ref. |35| . 

We begin our discussion of intra-ensemble decays by discussing the partial-width contribution arising from the 
decays of a given ax into multiple, lighter axion modes. The leading contribution to this partial width comes from 
the quartic interaction term appearing in Eq. (|2.26l) . which gives rise to decays of the form ax — >■ axiax20-\a, where 
axi, with i = {1, 2, 3}, are lighter axion states whose masses satisfy the constraint A > Ai + A2 + A3. The Feynman 
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rule for the corresponding four-point interaction vertex is 
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from which the contribution to V\ from this three-body decay is found to be 
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where now 
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(pi +P2)^ and TO23 = (p2 +P3)^7 with "pi representing the four-momentum of the final-state oa;- 



The limits of integration for the dm^2 integral are (?Tii2)max = (A — As)^ and (mi2)^j„ ~ (\\ + ^■if' , and since we 
are primarily interested in the parametric dependence of the partial width on A, Mr^^ etc., it will be sufficient for 
our purposes to construct an upper bound on T{ax -^ aAiaA2'3.A3) by setting ('Ti23)max ^^ (A — Ai)^ and (rn23)f^i^ ~^ 
(A2 -I- Aa)^. Doing so, we obtain the result 

r(aA^aA,aA.aA3) < ^ (|^) (A^^a)^ fll A^^a.) 

X TAi + A2 + A3 - X\ Ta^ - A2 + A^ - A^ + 2AA2 + 2X1X3] . (3.10) 

Note that the asymmetry of this expression under permutations of the Xi is due to the asymmetric limits of integration 
we have adopted in order to construct this bound. In order to obtain the total contribution Fa (a — > 3a) to the partial 
width of a given ax from decays of this form, we need to sum over the contributions from all kinematically allowed 
decays of the form a\ — >■ aAiaA2 0A3- Next, we approximate the sums over the different allowed final-state axions with 
integrals over dXi, dX2, and dX^. Furthermore, Eq. (|2.24p implies that X^A\ < \/2 for all A. Therefore, in order to 
obtain an upper bound on Fa (a — !• 3a), we make the replacements X^Ax^ — >■ V^ and X^Ax — > \/2. Doing so, we obtain 
our final result: 



rA(a-^ 3a) < 






(3.11) 



Given this result, we are now prepared to address the question of whether axion decays to other bulk axions can 
ever contribute significantly to the total width Fa of any aA . For example, the bound in Eq. p. lip implies that the 
ratio of Fa (a — >■ 3a) to the decay rate Fa (a — >■ 77) given in Eq. (JSHJ) is bounded from above by 
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As we shall see in Sect. IVH for 0(1) values of ^ and go, the phenomenologically preferred ranges for fx and Kq 
turn out to be fx ~ 10^^ - lO^^ GeV and Ac ^ 10^ - 10^ GeV, while M^ is bounded from below by Eq. (jXi)) . 
Within this parameter-space regime, we find that the ratio in Eq. p.l2p will be vanishingly small, as desired, unless 
A ^ 10^*^ GeV. Since this is far larger than the cutoff scale fx in this same regime, we conclude that decays of the 
form aA — > aAiaAa^As will not play a significant role in the phenomenology of realistic bulk axion models of dynamical 
dark matter. Thus such decays can be safely neglected in computing the total decay width of a given aA. 

Not only can the oa decay to final states comprising lighter axion modes alone, but they can also decay into final 
states which include other bulk states. In the minimal bulk-axion theory under discussion here, these include KK 
graviscalars and KK gravitons. (Note that the vector degrees of freedom /ij^g with n > in the gravity multiplet 



(0) 



do not couple to the aA in the linearized-gravity limit in the unitary gauge, and the zero-mode n\^ vanishes due to 
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the orbifold projection.) Therefore, we must also assess whether decay channels involving these KK gravitons and 
KK graviscalars can provide an appreciable contribution to F^^ '. In the five-dimensional theory under discussion 
here, in the unitary gauge, the only physical graviscalar present is a single radion mode, which we assume here to be 
sufhciently massive {e.g., as the result of some stabilization mechanism) as not to be relevant for ax decays. As for 
decays involving KK gravitons in the final state, a rough upper bound on their contribution to T\ for the case of a 
single, flat extra dimension will be given in Ref. |35| within the framework of linearized gravity. As wc will see, the 
leading contribution comes from two-body decays of the form a\ -^ GjTJax', where GfZ denotes the KK graviton 
with KK mode number n. The total contribution from decays of this sort, summed over all kinematically accessible 
combinations of n and A' , is found to be approximately [35| 
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where Kixx) and E{xx) respectively denote the complete elliptic integrals of the first and second kind, with xx = 
(A-A')V(A + A')2. 

In order to compare T{ax —J- Ga) to the rate for ax decays to SM fields, it is necessary to integrate Eq. p.l3p 
numerically, as a function of fx, Ag, and Mc. The results of such an analysis are given in Ref. [35|- Here, however, to 
illustrate our point, we simply choose a set of benchmark values typical of a phenomenologically consistent scenario 
for which r{ax —> Ga)/T{ax —>■ 77) is roughly maximal. Specifically, we take Ac = 1 TeV and Mc = 10~^^ GeV, with 
go = ^ = i. We then find that T(ax — ^ Ga) always remains several orders of magnitude smaller than T{ax — >■ 77) for 
all values of fx ^ 10^^ GeV. We therefore conclude that the decays of ax to KK gravitons will not have a significant 
impact on the total widths of the ax- 
Taken together, these results strongly suggest that intra-ensemble decays do not play a significant role in the 
phenomenology of bulk-axion dynamical dark matter, and can therefore be neglected. Therefore, from this point 
forward we will ignore intra-ensemble decays and identify Tx, as given in Eq. p.8p . with the total width of any given 
ax- 



Axion Lifetimes Across the Ensemble 



In Fig. [1] we show how the lifetime tx = 1/Fa of an axion mass eigenstate ax behaves as a function of A. The 
left panel shows the results for an axion with Aq = 1 GeV, while the right panel shows the results for an axion 
with Aq = 1 TeV. In each case, we have taken go = £, — ^ and set Mc = 10~^^ GeV. In each of the two panels, 
the three solid curves correspond to three different choices of fx for a photonic axion with c^ = 1. The solid 
red curve corresponds to fx = 10* GeV, the solid green curve corresponds to fx = 10^^ GeV, and the solid blue 
curve corresponds to fx = 10^^ GeV. The dashed curves correspond to the same choices of fx for an axion with 



1 and Ci = 1 for i = {e 



to the thresholds at 



■'ej i^lj.-, 



Vt}- The series of kinks which are evident in each dashed curve correspond 
and rUf^ above which new decay channels for ax open up. The sharp drop in tx at 
around A ~ 400 MeV is the result of the a -> 7r+7r~7r° decay channel opening up. 

One significant property of the decay rates of the ax in bulk-axion scenarios can be readily appreciated upon 
comparing the curves appearing in the two panels of Fig. [TJ the total width F> is independent of mx (and therefore 
independent of Aq) in the limit in which A ^ tttjiJ^/Mc- This can also be seen from Eq. (|3.8I) . It therefore follows that 
the corresponding curves appearing in the two panels of this figure should coincide for values of A above the threshold 
at which this condition is met for both of the selected values of A^. Indeed, we see that this is in fact the case. 
The fx = 10^^ GeV curves coincide for nearly the entirety of the range of A shown, since TTm\/Mc is approximately 
9.95 X 10-^2 Qev for Ac = 1 TeV and is far smaller for Ac = 1 GeV. The fx = lO^^ GeV curves, on the other 
hand, begin to coincide only for A ~ lO^'^ GeV, which is just above the threshold ■Kra\/Mc ~ 9.95 x 10~'* GeV for 
Aq = 1 TeV. 

There is, however, an even more important message to be gleaned from Fig. [TJ Note that in each panel, we have 
included for reference a set of horizontal, dashed lines indicating the time scales associated with the beginning of BBN 
(^BBN ~ 1 s), the present age of the universe (inow ^ 4.3 x 10^"^ s), and the usual limit on the lifetime r^ of a single 
decaying dark-matter candidate x given in Eq. (|1.2p . These benchmark times are absolutely critical for the survival of 
our dynamical dark-matter model- Any ax with a lifetime that falls between ieBN and r^ has the potential to disrupt 
BBN predictions for the abundances of light elements, distort the CMB to an unacceptable degree [2l, [Sa], produce 
too large a flux of X-ray or gamma-ray photons, etc- For this reason, the success of our dynamical dark-matter model 
rests upon the assumption that such ax have sufficiently small relic abundances Vix that these decays are harmless. 
It is in this manner that lifetimes must he balanced against abundances across our dark- axion towers- 
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FIG. 1: The lifetime tx of the axion mass eigenstate ax, shown as a function of its mass A for Ag ~ 1 GeV (left panel) and 
an axion with Ag = 1 TeV (right panel). In both panels, we have set ^ = ^g = 1 and have chosen Mc = 10~^^ GeV — 
a value just above the lower bound imposed by modified-gravity experiments. The three solid curves in each panel shown 



correspond to different choices of /x for a "photonic" axion with the coupling-coefficient assignments Cy 



1, Co 



0, and 



Ci = for all fermions i — {e, Ve, i^fj., Vt}- The solid red curve corresponds to fx ~ 10 GeV, the solid green curve corresponds 

to fx ~ W^^ GeV, and the solid blue curve corresponds to fx 

choices of fx for an axion with the coupling assignments c-y — Cg 

opening up of new decay channels as A is increased past a series of kinematic thresholds associated with decays to neutrino 

pairs, electron pairs, muon pairs, and ■k'^-k~-k^ . The horizontal lines indicate the time scales associated with the onset of BBN 



■ IQ^^ GeV. The dashed curves correspond to the same 
1 and Ci = 1 for all i. The kinks in the curves reflect the 



(teBN ~ 1 s) , the present age of the universe (tno 
T-^ of a single-particle dark-matter candidate. 



4.3 X 10 s), and the usual lower limit given in Eq. (|1.2|l on the lifetime 



The results in Fig. [T] also highlight another important aspect of our dynamical dark-matter ensemble, which is that 
at a given time i, only a fraction of the a\ — those which have not already decayed — can contribute significantly to 
^CDM- Therefore, since T\ increases monotonically as a function of A, there exists a maximum value Adoc for which 
ax may be considered stable for particular time scale t (in the sense that Txt < 1) and which is potentially capable of 
contributing significantly to J7cdm- For a photonic axion, for example, the approximate form of Adoc can readily be 
obtained in both the y <C 1 and y ^ regimes by inverting Eq. (|3.1[) : 
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In Fig. [2j we illustrate the behavior of Adec (solid curves) for a photonic axion a function of jx for several choices 
of M^. In each case, we have fixed Ag = 10 GeV and set ^ = (7g = c-y = 1. For reference, we have also included the 
corresponding curves for two other critical values of A in any given axion KK tower for each M^. These are the mass 
Ao (dashed curves) of the Hghtest axion mass eigenstate a\^ and the mass Atrans (dotted curves) defined in Eq. ()2.28p 
which dehneates the transition point between the small-A and large-A regimes. We have also included a pair of black 
curves indicating the asymptotic behavior of Aq and Adcc for large M^- Note that in this limit, mixing is negligible, 
and Ao ~ mx for all values of fx ■ 

The behavior of Aq, Adoc, and Atrans depends primarily on the value of the mixing parameter y. When fx is 
sufficiently large that y ^ 1 for a given value of Mc {i.e., in the lower right portion of Fig. [2]), we find that Aq ~ mx, 
as expected. At the same time, since mixing is negligible in this regime, increasing fx results in a uniform suppression 
of the decay widths of all ax, and hence Adoc increases with increasing fx (i.e., in the upper right portion of the 
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FIG. 2: Curves showing a variety of critical values of A as functions of fx for Ag = 10 GeV and several different choices of Mc 
ranging from 10~^® GeV to 10~* GeV. Each dashed curve indicates the mass Ao of the of the lightest axion mass eigenstate 
axQ for a given choice of Ma- Each solid curve indicates the mass Adec of the heaviest oa for that choice of Ma which has not 
decayed by present time, assuming a photonic axion with c^ = 1. Each dotted curve marks the transition point Atrans between 
the small- A and large- A regimes for the same choice of Ma, as defined in in Eq. 1)2. 28p . The black curves indicate the asymptotic 
behavior of Ao and Adcc for large Ma- 



figure). However, as fx decreases past the point at vi^hich the Aq and Atrans curves intersect, we pass from the j/ S> 1 
to the y ^ 1 regime. In this latter regime, mixing is significant and Aq ~ A/c/2. At first, Adcc still decreases with 
decreasing fx in this regime, as those modes with A < nm'j^/Mc continue to destabihze down the tower. However, 
as fx decreases still further, to the point at which the Adcc and Atrans curves intersect, the a\ in the A < nm'j^/Mc 
regime begin to destabilize as well. The dependence of T\ on fx is qualitatively different for such modes, as indicated 
in Eq. (|3.14[) . and consequently Adcc actually begins to increase with decreasing fx- These observations will turn out 
to be critical in interpreting the results to be derived in Sect. |Vl 



IV. CHARACTERIZING THE CONSTITUENTS: RELIC ABUNDANCES 

In the previous section, we focused on one of the two crucial properties of the particles that constitute a dynamical 
dark-matter ensemble: their stability. In particular, we examined the decay widths of the various a\ and investigated 
how these widths scale with A. In this section, we focus on the other property: their abundances. We begin by 
establishing a consistent cosmological context in which to situate our bulk-axion theory. We then proceed to evaluate 
the various mechanisms through which a population of a\ can be produced in the early universe, and explain why 
misalignment production is favored from the perspective of dynamical dark matter. This thereby justifies the emphasis 
placed on this mechanism in Ref. Jj. We then derive explicit formulae for fix, and demonstrate that the proper balance 
between T\ and Qx is indeed realized in the context of misalignment production. In the following section, we will then 
use these results for Tx and fi^ to characterize the aggregate properties of the entire ensemble, such as its equation 
of state, its total relic abundance, and the way in which that abundance is partitioned among its constituents. 



A. Standard and Low- Temperature Reheating (LTR) Cosmologies 



Before embarking on a discussion of axion production in the early universe, however, we must first specify the 
cosmology in which that production occurs. This is particularly relevant in the context of theories involving large 
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extra dimensions, since the properties of the early universe in such theories can differ dramaticahy from those which 
characterize the standard cosmology. For example, the pres ence of substantial energy density in the bulk can alter 
the expansion rate of the universe in a significant way [33], and late decays of KK excitations of the graviton (or 
other bulk fields) can disrupt BBN, produce visible distortions in the diffuse photon spectrum, etc. For this reason, 
such scenarios must obey stringent constraints [28| on the so-called "normalcy temperature" T*: the temperature 
below which the universe is effectively four-dimensional, in the sense that the bulk is essentially empty of energy 
density and the radii of the extra dimensions can be regarded as fixed. These bounds come from a diverse array of 
considerations and leave a very narrow window of 4 MeV ^ T* < 30 MeV for this normalcy temperature. The most 
attractive solution for arranging such a value of T* is to posit a very late period of cosmic inflation, precipitated by 
a brane- localized inflaton [28|, with a reheating temperature Trh — T^- 

If Tyih is indeed to be identified with T*, then the universe must be described not by the standard cosmology, 
but by an alternative framework commonly dubbed the low-temperature reheating (LTR) cosmology [3a |. In the 
LTR framework, inflation occurs very late, and the energy density of the universe remains dominated by coherent 
oscillations of the inflaton field (j). Such oscillations behave like massive matter down to very low temperatures — 
potentially as low as a few MeV. Consequently, the universe undergoes an additional epoch of matter domination 
(MD) at early times, which ends only once decays of cj) into SM fields in the radiation bath have depleted the energy 
density p^ to the point that p^ = Prad- This condition defines the reheating temperature Trh, which is determined 
solely by the decay width F^ of the infiaton: 



TjiH — 



90 



1/4 



, .^ . yTV^> (4.1) 

where g*(T) denotes the effective number of massless, interacting degrees of freedom at temperature T. At temper- 
atures T < Trh the universe is effectively radiation dominated (RD), and maps onto the standard cosmology. This 
modification implies that the relationship between time and temperature in the LTR cosmology is quite different from 
that obtained in the standard cosmology at early times. In particular, one finds that in any universe which underwent 
a period of cosmic inflation immediately prior to the RD era, the relationship between time and temperature is given 
by 
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where Tmre ^ ©(eV) is the temperature associated with the usual matter/radiation transition, at which the energy 
density of the universe once again becomes dominated by the contributions from dark and baryonic matter. (For our 
purposes, it will be sufficient to approximate the present, ACDM universe as matter-dominated.) In the standard 
cosmology, Trh is high enough that the universe will be radiation-dominated at all time scales relevant to axion 
dynamics, all the way down to the time scale ^mre associated with this transition. By contrast, in a LTR cosmology 
with Trh ~ O(MeV), much of the relevant dynamics will occur while the universe is still dominated by coherent 
oscillations of the inflaton field. Such a cosmological modification can have profound effects on axion dynamics, even 
in the case of a four-dimensional axion [39J, |40[ . 

Since our model necessarily involves large extra dimensions, constraints on T^ would seem to require that the 
universe be described by the LTR cosmology, rather than the standard cosmology, at early times. Indeed, for this 
reason, we shall adopt such an LTR cosmology in what follows. This will actually turn out to be advantageous, as 
the relic abundance of a light scalar generated via misalignment production can differ substantially between the two 
cosmologies [40|. To better facilitate comparison between the results obtained in these two cosmologies, we therefore 
find it instructive to present results for each in parallel. The reader should keep in mind, however, that the results 
obtained for an LTR cosmology with 4 MeV ^ T* < 30 MeV should be taken, in some sense, as the "true" ones, given 
the constraints on T*. 

In this paper, we shall be interested in values of the confinement scale Kq which range between roughly 10 MeV 
and 100 TeV. When operating in the standard cosmology, we shall assume that Ag <C Trh, so that confinement takes 
place within the RD era. Conversely, when operating within the LTR cosmology, we will assume that Kq > Trh, so 
that confinement takes place before or during reheating, when the universe is dominated by coherent oscillations of 
the inflaton fleld. Note, however, that for smaller confinement scales Aq < Trh, the results for the LTR cosmology 
are identical to those for the standard cosmology. 
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B. Axion Production Mechanisms 

Having now set the cosmological tableaux, let us begin our discussion of axion production in the early universe. A 
number of production mechanisms can, in principle, contribute significantly to the axion relic density. In generic 
models, three such mechanisms typically tend to provide the dominant contribution to the relic density of any 
individual axion field. One of these mechanisms is thermal production via the interactions of the a\ with the SM fields 
in the radiation bath. The other two are non-thermal in nature and can generate a population of cold axions. These 
are misalignment production and production from the decays of topological defects (in particular, cosmic strings) 
associated with the breaking of the global U{l)x symmetry. 

We begin with a discussion of thermal production. A number of processes can contribute appreciably to thermal 
axion production in the early universe, depending of course on the magnitudes of the couplings between the axion 
in question couples to the SM particles. Among hadronic processes, qg — ?> qax, qq -> ga\, gg —?> ga\, etc., dominate 
for T > Aqcd, while pion-axion conversion off nuclei (including all processes of the form Nn — )■ N'ax, where 
N,N' = {n,p}) and the purely pionic process tttt — ^ nax dominate at lower temperatures. Since we are interested in 
values of M^ which are far below Trh, it follows that T >> A for a large number of the a\ for at least some of the 
post-inflationary era. These a\ can therefore be considered effectively massless at such temperatures. In this massless 
limit, the rate for each of the axion-production processes enumerated above (except for the inverse-decay process, 
which is generally sub-leading) takes the rough, parametric form 



^'-(A^Aa)2 . (4.3) 
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In other words, since kinematic distinctions between states for which T 3> A are unimportant, the dependence of their 
production rates on A occurs primarily through the coupling factor X'^A\. This factor too is effectively independent 
of A when A > ttttiJ^/Mc', hence, at a given temperature T, any a\ with a mass A in the range T:m\ /Mc ^ A < T will 
be produced at essentially the same rate. Moreover, production rates actually increase with increasing A for those 
modes with masses A < Trmj^/Mc- This means that at a given temperature T, the heavier ax in this mass range will 
actually be produced from the thermal bath at an equal or higher rate than the lighter a\ — at least until A becomes 
comparable with T, and their production rates are Boltzmann suppressed. Since Tx also increases with A, the less 
stable states will be thermally produced at an equal or higher rate than the more stable ones. Such a relationship 
is clearly undesirable in models of dynamical dark matter. Moreover, the majority of light axions produced through 
interactions with the thermal bath would be relativistic at the time of production, and therefore not cold. 

From these considerations, we conclude that if an ensemble of axion KK modes is to constitute the majority of the 
dark-matter relic density, thermal production must contribute only a negligible fraction of the total relic abundance 
of each ax, with the remainder of that abundance generated through non-thermal means. This requirement places 
a non-trivial constraint on scenarios of this sort, a detailed analysis of which appears in Ref. [35|. We will defer 
the discussion of how this constraint restricts the parameter space of our model until Sect. IVII For the moment, we 
simply note that this constraint exists, and proceed to discuss non-thermal mechanisms for axion production. Note, 
however, that in traditional models of KK dark-matter (either single-component [7] or multi-component [6]), in which 
the dark-matter candidates are stable, thermal freeze-out can be a viable production method for generating relic 
abundances. 

One method in which a non-thermal population of axions may be generated in the early universe is production 
via the decay of cosmic strings associated with the broken global U{l)x symmetry. However, this mechanism can 
contribute significantly to axion production only if Hj > fx , where Hi is the value of the Hubble parameter during 
inflation, so that those cosmic strings are not inflated away. Since the value of Hi is relatively unconstrained, and 
since astrophysical and cosmological constraints will turn out to require fx to be quite large, in the rest of this paper, 
we shall assume that Hi <C fx- We will therefore not consider axion production from U{l)x string decay. However, 
it should be noted that in other scenarios (or in other regions of parameter space) , axions produced via cosmic-string 
decay could have important phenomenological consequences in dynamical dark-matter models, and this production 
mechanism therefore deserves further study. 

Finally, we turn to non-thermal axion production via the misalignment mechanism. As we shall see, this turns out to 
be the most promising axion-production mechanism from the perspective of dynamical dark matter. The basis of this 
mechanism is that at temperatures T ^ Aq , the instanton- induced contribution to the axion potential in Eq. (j2.18p 
effectively vanishes. This implies that the only contribution to the axion mass matrix at such high temperatures are 
the diagonal contributions from the KK masses: no mixing occurs, and consequently the mass eigenstates are merely 
the KK eigenstates a„. While the potential for each a„ with n 7^ is therefore no n- vanishing, due to the presence of 
the KK masses, and is minimized at a„ = 0, the potential for the zero mode oq vanishes. In the absence of a potential 
for ao, there is no preferred vacuum expectation value (ao) which minimizes V{ao); indeed, any (ag) < fx is as good 
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as any other. This means that when the U{l)x symmetry is broken, the value of (ao) within a given domain is 
essentiahy arbitrary. Thus, immediately after this phase transition occurs, one would expect to find a set of domains, 
each with a different homogeneous background value for the axion field, which would generically be expected to be 
^(/x)i but could in principle be smaller. Our ignorance of this initial value of (oq) is commonly parameterized by 
a "misalignment angle" 9, so that the initial conditions at the time at which the U{l)x symmetry is broken can be 
written as J2Q | 

(ao) = dfx , (a„) =0 for n ^ . (4.4) 

Indeed, as was noted in Ref. [li], this initial condition follows from U{l)x invariance, which manifests itself here in 
the form of a five-dimensional shift symmetry under which a ^ a + c, where c is a constant. 

At lower temperatures, however, the situation changes, as instanton effects generate a brane mass mx{T) for the 
bulk axion. Here, we write mxiT) rather than mx in order to emphasize that this instanton-induced mass term is 
temperature dependent, and reserve the symbol mx (without the argument) to refer to the constant, late-time (i.e., 
low-temperature) value of mx{T). Assuming that the instantons associated with the group G behave analogously to 
QCD instantons, Eq. ()2.5|) implies that mxiT) scales roughly like (Ag/T)*. Thus, when T ^ Kg, the off-diagonal 
terms in the mass matrix A^^,„j(r) become appreciable, and the a„ are no longer mass eigenstates. In this regime, 
the equations for a„ form a coupled system [26| , with the evolution of each such field governed by an equation of the 
form 

oo 

a„ + 3i/d„ + ^rA(r)[/A„(r)a„ + ^ Ml„,{T)a,n = , (4.5) 

A m=0 

where U\n{T) denotes the unitary matrix in Eq. (j2.23p . a dot denotes a derivative with respect to the time t, and 
H is the Hubble parameter. Note that since the mass eigenvalues A(r), decay widths rA(r), and even the rotation 
matrix U\n{T) itself all depend on m,x, these quantities all implicitly depend on temperature, and hence on t. 

During any period in which Tx{T) and A(r) vary appreciably in time, it is in general not possible to write down 
an exact, closed-form solution to the coupled system in Eq. (j4.5|) . Fortunately, however, m,x{T) can be regarded as 
effectively constant during most of the history of the universe. The only exception occurs at temperatures around 
T ~ Ag, where mx{T) rapidly rises from a negligible initial value to the asymptotic value it attains at T ^ Kg- At 
all other times, m,x{T) is well approximated either by zero or by tox, and the system of equations therefore decouples 
in the mass-eigenstate basis. In this basis, the time-evolution of each field a\ is governed by an equation of the form 

ax + ^flA + Taoa + A^flA ^ . (4.6) 

In arriving at this expression, we have used the fact that within an RD or MD era, H is approximately given by the 
relation 3iJ ~ w/i, where 

[3/2 iuRD ,^ ^, 

K = { ' (4.7) 

]2 iuMD. ^ ' 

Exact, closed-form solutions for a\ and a\, given an evolution equation of this form, do exist, and we present these 
solutions in the Appendix. It may be observed, however, that Eq. (14.61) is simply the equation of motion for a damped 
harmonic oscillator with a time-dependent damping term. As we shall see at the end of this section, it turns out 
that Fa ^ 3-ff at the time when A ~ 'iH/2 for all ax. It then follows that the solutions for each ax can be divided 
into two regimes, depending on the relationship between A and H at any given time t. When A < 37J/2, ax does 
not oscillate, and therefore its energy density scales approximately like vacuum energy. By contrast, when A > ?)H /2, 
ax oscillates coherently around the minimum of its potential, with oscillations damped by a "friction" term with 
coefBcient (3iJ + FA). 

Despite the fact that m.x{T) is effectively constant both well before and well after the time scale Ig at which T — Kg, 
the non-trivial dynamics of the ax ait ^ Ig certainly can play a crucial role in establishing the initial relic abundances 
for these fields. Nevertheless, while such a time-dependence can indeed have a significant quantitative impact on the 
results of relic-density calculations in certain cases, as has been shown to be the case with a standard, four-dimensional 
axion [22] , a great deal of information can be obtained by working in the "rapid-turn-on" approximation, in which 
we approximate 

mx{t) = mxQ{t-tG) . (4.8) 

In this approximation, the procedure for calculating the background value of each ax at any time t is clear [l|. At 
early times, when t < Ig, mx{T) = 0, and hence the a„ remain fixed at the initial values given in Eq. (14.41) . At 
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t = tc, rnx{T) immediately assumes its constant, non-zero, late-time value m^, and each of the mass eigenstates a\ 
acquires a background value proportional to its overlap with oq: 

(oa) = OfxAx , {ax) =0 at t ^ to ■ (4.9) 

Even though all of the a\ acquire background values at t = to, only those fields for which A > 'iH{tc)/'2, begin 
oscillating immediately at the time of this phase transition. As discussed above, all other, lighter a\ will begin 
oscillating later, once the A > 3iJ/2 threshold is crossed. The time tx at which a given a\ begins to oscillate is 
therefore given by 

tx = max{^,iG} , (4.10) 

where kx is the value of k corresponding to the epoch during which this oscillation begins. At times to < t < tx, a 
given ax continues to behave like vacuum energy rather than like matter, and thus properly contributes not to ricDM, 
but to the dark-energy abundance. (Note that this definition of tx is slightly different from the one given in Ref. [l| , 
where tx was simply defined as kx/X, regardless of its relationship to to.) 

For any given ax, however, the relevant quantity for dark-matter phenomenology is not the value of ax itself, but 
its energy density px, which is related to ax and dx by the relation 

Px = IH + >^^4] ■ (4.11) 

At early times, when t < to, only oq has a non-zero background value in the rapid-turn-on approximation, and since 
this field is massless, its energy density vanishes. At t ^ to, however, Eq. (|4.9p implies that each field acquires an 
initial energy density 

Pxitc) = lo'f'x\'Al . (4.12) 

Since ax remains effectively constant until t — tx for any field for which tx > tQ, we also see that px{tx) ~ Pa(^g)- 
This implies that the energy density stored in such a field behaves like vacuum energy until t = tx, at which point 
the field begins to oscillate coherently around the minimum of its potential. The energy density stored in such 
oscillations, as is well known, scales like massive matter. At late times t ^ tx, when the time scale associated with 
these oscillations becomes rapid compared to the time scale over which the amplitude of ax changes appreciably, the 
virial approximation implies that px ~ (^a)' "^h^re (of) denotes the average of d^ over one cycle of oscillation. In 
this regime, one finds that 

Px{t) - Pxita) (jX^-^'^'-'''^ (4.13) 

for each px during the epoch in which oscillation began, with pxita) given in Eq. ()4.12|) . Computing px during 
subsequent epochs is simply a matter of applying Eq. (|4.13l) iteratively with the appropriate boundary conditions for 
Px at the transition points at which k changes. 

All that remains, then, in order to specify the energy density px associated with a given ax for any particular choice 
of model parameters is to determine the time scales tQ and tx as a function of those parameters. Indeed, the results 
for Px clearly depend sensitively both on the epoch during which abundances are established, and the epoch during 
which oscillation begins. In principle, tQ could fall within the RD era, during the usual MD era, or during reheating; 
likewise, tx could occur at any time at or after tQ. However, as will be shown in Sect. |Vl Kq ^ 10 MeV is required 
in order for our ensemble of ax to yield a relic abundance on the order of ricDM- For such values of Kq, tQ ^ tMRE. 
Moreover, for 0(1) values of qq and ^, tx <C tMRE as well. Therefore, when discussing the standard cosmology, we 
will focus on the case in which tQ < tx < tMRE for all ax- Furthermore, when operating within the context of the 
LTR cosmology, we will implicitly assume that tQ < tx < tRH for all ax, so that all fields begin oscillating while 
the energy density of the universe is still dominated by coherent oscillations of the inflaton field. This is justified by 
the fact that tRH ^ 10^^ — 10^* s for reheating temperatures Trh within the phenomenologically allowed window 
4 MeV ^ T^, < 30 MeV. In conjunction with the experimental bound on Mc given in Eq. (|3.4[) . this implies that 
tA ^ tRH for all ax in any given tower, unless Aq < 10 MeV. In summary, we shall therefore assume that 

standard cosmology : tRH < to ^ tx < tMRE 

LTR cosmology : tp <tQ <tx < t^H (4-14) 

in what follows, where tp is the Planck time. Note that in the standard cosmology, tj^H is assumed to be so early 
that all tQ of interest will easily satisfy the lower bound. Also note that in the LTR cosmology, modes for which tx 
occurs during or prior to inflation will infiate away and therefore carry zero abundance at present time. 
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Axion Relic Abundances 



We now provide explicit expressions for p\ in both the standard and LTR cosmologies. We begin by considering 
the case of the standard cosmology, in which k — 3/2 at all relevant time scales prior to matter-radiation equality, 
and K = 2 after the transition to matter-domination at ^mre- In this cosmological framework, given the regimes for 
to and t\ specified in Eq. (|4.14p . it therefore follows that 
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(4.15) 



By contrast, in the LTR cosmology, the energy of the universe remains dominated by coherent oscillations of the 
inflaton field from the end of inflation until a very late time ^rh ^ 1 s. Thus, if the axion fields begin oscillating at 
a time t\ < ^rh, as specified in Eq. (I4.14p . we initially have k — 2, followed by a transition at tjm to the usual RD 
era, in which k = 3/2. This signifies that in the LTR cosmology, we have 
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(4.16) 



In other words, Eq. (|4.16l) replaces Eq. (|4.15l) in the context of the LTR cosmology, in which the usual RD era is 
preceded by an initial period of matter domination. It is worth emphasizing here that the value of tx for a given mass 
eigenvalue A will, in general, differ between the two cosmologies, due to the differing relationship between H and t at 
times t < fRH- 

Comparing Eqs. (|4.15l) and (|4.16p . we see that the cosmological context in which the axion fields evolve can have a 
potentially dramatic effect on the late-time results for the various px- However, for any given mode, Eq. (|4.10|) implies 
that the magnitude Sltr of that suppression depends on the relationship between A and to- Comparing Eqs. ()4.15p 
and (|4.16l) . we find that the axion energy densities are suppressed in the LTR cosmology, relative to the standard 
cosmology, by a factor 
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These results imply that the energy-density contributions from those modes which begin oscillating at Iq are sup- 
pressed in the LTR cosmology, relative to their value in the standard cosmology, to a greater degree than the contri- 
butions from those modes which begin oscillating later. These results are analogous to those obtained in Ref. [40|] for 
a standard, four-dimensional axion, where the suppression factor is referred to as yLTRy^std^ 

It is also possible (and indeed when Kq is large, more or less inevitable) that in the LTR cosmology, a great 
many of the heavier ax will begin to oscillate either during or prior to the end of inflation. Since the scale factor R 
grows exponentially during this epoch, the energy density in stored any such mode, which scales like px oc R^^ ^ will 
effectively be diluted into irrelevance by this rapid expansion. As long as Ag < Tmax, *-e-, as long as confinement 
occurs only after inflation is over, no energy density is stored in these modes during inflation. As a result, the energy 
density is given by Eq. (I4.16|) as usual. However, if conflnement occurs before or during inflation, px for any mode for 
which A > 3Hj/2 will be exponentially damped by Hubble dilution, and it is therefore reasonable to take px — for 
any such mode. 

From the results in Eqs. (|4.15l) and (|4.16p . it is straightforward to obtain the relic abundance fix = px/Pait for 
each ax- Let us begin by addressing those modes for which tx — to- Since the critical density for a flat universe is 
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given by pcrit = SH'^Mp, we find that in the rapid-turn-on approximation, the contribution to the dark-matter rehc 
abundance from each such mode at a given time t in the standard cosmology is 
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By contrast, in the LTR cosmology, the corresponding result is 
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For the rest of the ax {i.e., those for which tx > to), the corresponding results in the context of the standard cosmology 
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whereas in the context of the LTR cosmology, we instead have 
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Note that this is valid only for those ax for which A < 3Hj/2. For those modes with A > 3Hj/2, as discussed above, 
effectively fix = due to Hubble dilution during inflation. During periods of matter domination, we see that, px and 
Pcrit scale identically with time, and consequently fix remains constant in each of these expressions. During periods 
of radiation domination, on the other hand, px falls faster with time than pcrit, and fix grow like t^'^ . 

Note that for t = t„ow, Eqs. (|4.18p through (|4.2ip take the forms specified in Eqs. (54) through (60) in Ref. Q, 
where the brane mass m appearing in these equations is identified with mx for a bulk axion. In other words, the 
distinction between the fix expressions for the standard and LTR cosmologies given here is tantamount to identifying 
the era in which the initial abundances are established in Ref. [1|. 

Since we will be primarily interested in scenarios in which the dark-matter relic abundance ilcDM receives contri- 
butions from a large number of different axion mass eigenstates, it is of critical importance to determine precisely how 
fix scales with A, as this describes the relative contributions to the total relic abundance i^cDM from the various states 
in the tower. Fortunately, in the rapid-turn-on approximation, the dependence of fix on A stems from only three 
factors, as discussed in Ref. jl|. The first, of course, is the choice of cosmology. The second is how A compares to the 
mass scales mx and ttto^/Mc, which determines the overlap between ax and the KK zero mode ag, and therefore the 
initial displacement of a^ at i = tx. For those modes for which A > inax{?7ix, 7rm^/Mc}, the A^/m^ term dominates 
in the factor in brackets appearing in each of Eqs. (|4.15l) and (|4.16p . and fix acquires a factor of A~^. By contrast, for 
those modes for which A < maxJTTix, ttto^/Mc}, the constant terms dominate, and fix acquires no such dependence. 

The third factor which determines how fix scales with A in the rapid-turn-on approximation is whether or not the 
axion mode in question begins oscillating at tc, or at some later time. Indeed, in Ref. Jj, these cases were referred 
to respectively as the "instantaneous" and "staggered" turn-on regimes. Comparing the expressions valid for tx = to 
in Eqs. (|4.18l) and (I4.19P with those valid for tx > ta in Eqs. (|4.20p and (|4.2ip . we see that the fix for those ax which 
begin oscillating after the confining phase transition takes place acquire an additional dependence on A. The precise 
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form of this dependence depends on the cosmological context within which the model is embedded: in the standard 
cosmology, it is A~^/^; in the LTR cosmology, it is A~^. Physically, this factor stems from the fact that prior to the 
time it begins to oscillate coherently, the energy density in any given ax remains constant, whereas after oscillation 
begins, it scales like massive matter. Therefore, the later a given mode begins to oscillate, the longer the energy 
density stored in that mode will remain unaffected by cosmic expansion, and therefore the larger the present-day 
value of ft\ will be. 

In order to illustrate the implications of these effects, in each of the panels of Figs. [3] and S] we track the evolution 
of fix for a representative sample of ax within a given theory from to to present time. The curves shown in Fig. [3] 
reflect typical results which arise in the context of the standard cosmology. The left panel corresponds to a scenario 
with a small confinement scale Aq = 1 MeV, a moderate value fx = 10^ GeV for the effective four-dimensional 
C/(l)jf -breaking scale, and a small misalignment angle 6 sa 0.04. The right panel corresponds to the opposite case: a 
scenario in which Aq = 1 TeV, with /j^ sa 2 x 10^ GeV and 9 = 1. In each scenario, we have taken ^ = ^g = 1, and set 
the compactification scale to be Mc = 10"^^ GeV. The curves shown in each panel (from top to bottom) correspond 
to the lightest two values of A, here referred to as Aq and Ai, corresponding to that particular choice of fx and Mc 
(Ao « 6 X 10-1^ GeV and Ai « 10"" GeV for the left panel; Aq « 5 x lO^^^ QeV and Ai w 2 x 10"" GeV for the 
right panel), along with several additional values of A, including A = {10~^°, 10~^, 10~^, 10""', 10~^} GeV. It should 
be noted that in the right panel, all curves shown lie essentially on top of one another, and are thus not individually 
apparent. It is worth remarking that y 3> 1 in the scenario depicted in the left panel, while y ^ 1 in the scenario 
depicted in the right panel. Also shown are a horizontal, dot-dashed line indicating the value for f2cDM observed by 
the WMAP experiment Q, as quoted in Eq. (|l.ip . and a pair of vertical lines indicating the positions of ^rh and 

^MRE- 

Both of these scenarios yield a total dark- matter relic abundance which is consistent with the WMAP results given 
in Eq. (|l.ip . However, it should be noted that the parameter assignments used for these figures have been chosen 
exclusively for purposes of illustration. In Sect. IVI[ we will discuss the situation using the parameters which are 
consistent with all observational and phenomenological constraints. 

Note that the sets of curves shown in the two panels of Fig. [3] differ quite significantly. For example, in the left 
panel, the effect of the non-zero decay widths of the heavier ax is readily apparent. Indeed, the curves corresponding 
to masses A > 10~^ GeV rapidly drop to zero at a time scale t ^ tx- By contrast, in the right plot, which corresponds 
to a scenario with y ^ 1, the decay rate of each ax is suppressed by a factor of A'*^^, which can be quite small in 
such a scenario. Consequently, all of the ax for which fix curves are shown in the plot are stable on cosmological time 
scales. 

The abundance curves shown in Fig. 31 on the other hand, display typical results obtained in the LTR cosmology. 
More specifically, the results shown here correspond to the parameter assignments Trh = 5 MeV and Hj = 10 GeV. 
As in Fig. |31 the two panels shown in this figure correspond to two different choices for fx, Mc and Aq which 
both yield a total present-day relic abundance consistent with WMAP data. The curves shown in the left panel 
correspond to a scenario with fx = 10^ GeV, Mc = 4 x 10~^^ GeV, and Aq w 37 MeV. This scenario exemplifies 
the case in which fx and Aq are both comparatively small. By contrast, the curves displayed in the right panel 
correspond to a scenario with far larger values for these parameters: /x = 6 x 10^^ GeV and Aq = 1 TeV, with 
Mc = 10^^ GeV. In both cases, we have taken £^ = gc = d = I. Once again, the curves shown in each plot 
correspond to A = {10-i°, 10-^ 10-^ lO'^, IQ-^} GeV, as weU as Aq and Ai. For the left plot, Aq = 2 x lO'^^ GeV 
and Ai = 6 X IQ-^^ Q^y. foj. ^j^g right plot, Aq = 5 x lO'^^ GeV and Ai = 2 x 10"" GeV. Again, in the left panel, 
the Ao and Ai curves are not apparent because they lie directly beneath the A = 10"^° GeV curve. 

The contrasting features between the two panels in Fig. 0] are predominately due to the differences between their fx 
and Ag values. First, as stated above, fx is quite small in the scenario displayed in the left panel, and consequently, 
the couplings between the ax and the fields of the SM are quite large. This implies that the decay rates associated with 
the heavier ax will be sizable in this scenario, and that a large number of these heavier modes will decay before present 
time. Indeed, the precipitous drop in each of the curves corresponding to a mass eigenvalue in the range A < 10~^ GeV 
in this plot is a consequence of the decay of these modes to SM fields. By contrast, in the scenario displayed in the 
right plot, fx is large enough that all of the ax for which fix curves are shown are stable on cosmological time scales, 
and no such effect is apparent. Second, A^ is also quite small in the scenario displayed in the left panel, and the 
confinement time scale to ~ 10"^ s is consequently quite late. As a result, tx = ic for all ax- This implies not only 
that all of the modes begin oscillating at the same time, but moreover, that the Qx only become non-zero quite late. 
By contrast, in the right panel, Aq is large and to is correspondingly quite early. This results in a situation in which 
the ^A for the lighter modes are staggered in time. As discussed in Rcf. [1], the primary consequence of this staggering 
is that the fix curves for these lighter modes depend more sensitively on A. 

The most important implication of Eqs. (|4.18p through (I4.21[) . however, is that fix decreases with increasing A 
regardless of the details of the cosmological framework. On the other hand, we saw in Sect. |TTT] that Tx increases 
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FIG. 3: The individual relic abundances Q,\ associated with a variety of different a\ in the standard cosmology, shown as a 
function of time i, for two scenarios with different values of the parameters /x, Ag, and Q. The left panel corresponds to a 
choice of /x r; 10®, Ag = 1 MeV, and Q « 0.04, while the right panel corresponds to a choice of /x « 2 x 10^ GeV, Ag = 1 TeV, 
and Q = \. In both cases, we have taken ^ = (7g = 1, with M^ = 10"^'^ GeV. The range of t displayed in each panel spans 
from the corresponding confinement time scale to for the hidden-sector gauge group G, and a vertical, dashed line indicating 
the time scale associated with matter-radiation equality has also been included for reference. The horizontal, dash-dotted line 
indicates the total observed dark-matter relic abundance, as measured by the WMAP satellite. It should be emphasized that 
in both of these scenarios, the total present dark-matter relic abundance contribution from all of the a\ in the tower reproduces 
this observed value of rJcDM to within the limits quoted in Eq. (|l.ip . 
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FIG. 4: The individual relic abundances Q,\ associated with a variety of different a\ in the LTR cosmology, shown as a function 
of time i, for two scenarios with different values of /x, Mc, and Ag. The left panel corresponds to a choice of jx 



Mc = 4 X 10" GeV, and Ag ^ 37 MeV, while the right panel corresponds to a choice of /x f^ 6 x 10 GeV, Mc 



: 10°, 
10"" GeV 

and Ag = 1 TeV. In both cases, we have taken ^ = go = Q = \^ with Trh = 5 MeV and Hi = 10 GeV. The range of t 
displayed in each panel spans from the end of cosmic inflation to present day, and vertical (dashed) lines indicating the time 
scales associated with the end of reheating and with matter-radiation equality have also been included for reference. The 
horizontal, dash-dotted line indicates the total observed dark-matter relic abundance, as measured by the WMAP satellite. It 
should be emphasized that in both of these scenarios, the total present-day dark-matter relic abundance contribution from all 
of the a\ in the tower reproduces this observed value of SIcdm to within the limits quoted in Eq. (|l.ll) . 



24 

monotonically with A. This observation is indeed encouraging, in that it suggests that Tx and 17a possess the 
appropriate, reciprocal relationship needed for an ensemble of a\ to serve as dynamical dark matter. In the following 
section, we will quantify more precisely the relationship between T\ and fix and demonstrate that this is indeed the 
case. The results of the present section therefore attest that misalignment production, in stark contrast to thermal 
production, is an ideal mechanism for the generation of axion relic abundances in dynamical dark-matter models. 
Moreover, as we shall soon demonstrate, this mechanism dominates in the regime of model-parameter space in which 
an ensemble of ax tends to be phenomenologically viable, in the sense that it correctly reproduces the observed 
dark-matter relic abundance, while at the same time satisfying all relevant constraints from experiment, astrophysical 
observation, and cosmology. 

Before we proceed to analyze the collective properties of such ensembles, however, two brief comments are in order. 
The first of these concerns the validity of the rapid-turn-on approximation. As we have stated above, Eq. (14.61) is 
strictly valid only when A and Tx are essentially independent of temperature. However, there turn out to be certain 
situations in which the time-dependence of A and Tx at temperatures T > Aq is physically unimportant, and in which 
these quantities can be reliably treated as constants throughout the period in which ax are oscillating. One such 
situation arises in cases in which tx > to for all of the ax which contribute meaningfully to the total dark-matter 
relic abundance, and therefore coherent axion oscillations do not occur until after mx{T) is effectively constant. This 
situation tends to arise either when Aq is quite large, in which case mx{T) attains its constant, late-time value very 
early, or else when Aq is fairly small, but fx is quite large, in which case mx itself is extremely small. As we shall 
see in Sect. |Vl these turn out to be precisely the situations in which the total relic- abundance contribution from the 
ensemble of ax successfully reproduces the observed value for ficDM quoted in Eq. (jl.ip . This retroactively justifies 
our use of the rapid-turn-on approximation. 

Our second comment concerns the assumption that Tx is sufficiently small that the solution of Eq. (j4.6p for any 
given ax includes a period during which this equation of motion is effectively underdamped. This is critical, since 
in the absence of such a period, coherent oscillations cannot occur. Indeed, the energy density px stored in any ax 
for which Tx > 2A would never scale in an appropriate manner for that field to behave like massive matter; hence 
it would never contribute to r^cDM- However, it is not difficult to demonstrate that this situation essentially never 
arises in realistic bulk-axion scenarios, even for the most massive modes in a given tower. For example, consider the 
case of a purely photonic axion with c^ = 1. In this case, it follows from Eq. p.ip that the solution for ax becomes 
critically damped at a value Acd, which is determined by the condition 
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Solving this equation for Acd, we find that 
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which implies that Acd > /x/\/2G^. However, since the effective description of the theory in terms of a tower of 
axion modes breaks down at the five-dimensional C/(l)x-breaking scale fx <C fx, we are assured that Tx ^ 2A for 
all modes in such a tower. Indeed, this qualitative result is not specific to a photonic axion, but applies broadly to 
any axion field which couples to the SM fields with 0(1) coupling coefficients. As a corollary, this result also implies 
that the standard oscillation criterion A ~ 3H/2 will always be met before the decay criterion Tx ~ H. This indicates 
that indeed H, rather than Tx, sets the time scale at which oscillations begin, and that any given ax decays during a 
time frame in which its energy density can legitimately be described by Eq. (|4.13p . 



V. CHARACTERIZING THE ENSEMBLE: TOTAL ABUNDANCES, 
TOWER FRACTIONS, AND EQUATIONS OF STATE 



In the previous two sections, we derived expressions for the decay widths and relic abundances for the individual 
mass eigenstates ax in a mixed tower of KK axions. We have shown that these quantities scale with A in an appropriate, 
reciprocal manner for an ensemble of such states to serve as dynamical dark matter. We are therefore finally equipped 
to address the dark-matter phenomenology of the ensemble as a whole. 

As discussed in Ref. [ll, the crucial quantities which characterize a given dynamical dark- matter ensemble are the 
total relic abundance fltot, the tower fraction 77, and the effective equation-of-state parameter WeS- In this section, we 
investigate how these three quantities depend on the scales fx, Mc, and Aq which characterize a given bulk-axion 
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model and thereby assess which regions of parameter space are interesting from a dynamical dark-matter perspective. 
In the next section, we discuss the applicable phenomenological constraints on the model and demonstrate that 
substantial regions of parameter space exist within which all such constraints are satisfied. 

A. General Definitions 

The first of the three principal quantities mentioned above which characterize any given dynamical dark-matter 
ensemble is fltot- This is simply the total contribution to ficDM from all constituent modes in the ensemble which 
have already begun oscillating: 



l^tot = J^f^A . (5.1) 



The second is the so-called "tower fraction" t], which is a measure of how the total abundance J7tot is distributed 
across the ensemble. Specifically, 77 is defined for a given dynamical dark-matter ensemble to be the fraction of fltot 
provided by all of the oscillating components of that ensemble except for the one which yields the largest individual 
contribution. Explicitly, 

-1 "max /_ „N 

'7 = 1- 7^ — : (5-2) 

"tot 

where fimax = rnaxAJOA}- After all the a\ have begun oscillating, the lightest mass eigenstate oaq always yields 
the largest individual relic abundance V,\g, and therefore ^max = ^Aq- When 77 <C 1, essentially the entirety of fltot 
is provided by a single field, as in most traditional dark-matter models. By contrast, having 77 ~ 0{1) signals a 
departure from this traditional setup, which is indeed one of the hallmarks of dynamical dark matter. (Note that 
when we say that 77 should differ significantly from zero, we are willing to accept 77 ^ 0.1, as such values could result 
in observable differences from traditional models, but not, for example, 77 ~ 10"'^.) 

In our dynamical dark-matter framework, both fltot and 77 are intrinsically dynamical quantities, with non-trivial 
time dependences. For this reason, we will designate their present-day values as fi^Qt = ritot(^now) and 77* = 77(i„ow) 
in what follows. 

As discussed in Ref. [1], a given dynamical dark-matter ensemble as a whole can also be described in terms of a 
single, effective equation-of-state parameter Wcs: 

_ (I dlnptot , -,\ ,^ „^ 

where ptot = ^totPcrit- Indeed, one of the hallmarks of this framework is that even WaS itself is continually changing 
in time. We shall therefore define w^ = 7i'off(inow)- As discussed in Ref. [l|, this quantity is given by 

w* ~ n — n; (5.4) 

^^ 'tot '-now 

for any given dynamical dark-matter ensemble in which the widths, abundances, and densities of states obey the 
approximate scaling relations Vt w AF" and nr ~ BT^ , where nr denotes the density of states per unit decay width. 
Taken together, fitot , ??, and Wcs serve to characterize any given dynamical dark-matter ensemble. For the remainder 
of this section, then, our task will be to investigate how these three quantities depend on the parameters fx, Mc, and 
A(3 which characterize our bulk-axion model. Also recall that in Eq. (|2.2ip we defined the mixing parameter 

which quantifies the extent to which the different modes in the KK tower mix with each other. We will also therefore 
keep track of the corresponding values of y in our analysis. Moreover, as we have seen in Sect. IIV[ our results will 
also depend on the cosmological framework adopted. We shall therefore derive results in the context of the standard 
and LTR cosmologies independently. However, as discussed above, constraints on T* in theories with large, flat 
extra dimensions provide a strong motivation for working within the context of an LTR cosmology with a reheating 
temperature Trh ~ 0{MeV). For this reason, the results obtained for the LTR cosmology are more likely to be 
realistic. 
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Needless to say, phenomenological consistency imposes certain constraints on the parameters iltot, Vi ^^.d ""^off- For 
example, WMAP data require that fi^^j ~ ^cdm; likewise, w^ should not differ too significantly from zero. Beyond 
this, however, 77 and WefF are fairly unconstrained by data. Nevertheless, while any values for rj and Weff can be realized 
within the general dynamical dark-matter framework, we are particularly interested in situations in which -q^ is also 
significantly different from zero, for these are the situations in which our dynamical dark-matter ensemble represents 
a significant departure from traditional, single-component models of dark matter. 

B. Dark Towers: Relic Abundances and Tower Fractions 

In Fig. [5l we show how the total present-day dark-matter relic abundance ri^ot depends on fx, Mc, and Aq in the 
standard cosmology, assuming a photonic axion with c^ = 1 . Each panel in the figure displays contours of ^l^t ^^^ ^ 
different choice oi Aq. A dashed blue line highlighting the iljot — 1 contour has also been included in each panel. The 
red lines are contours of y: the solid red line corresponds to ?; = 1, which roughly indicates the transition point between 
the strongly-mixed regime (below and to the left of the contour) and the weakly-mixed regime (above and to the right 
of the contour). Proceeding from left to right, the dotted lines correspond to the values y = {0.01, 0.1, 10, 100}. In 
Fig. [51 we present the corresponding contour plots for the tower fraction 77*. Moreover, to complement the results 
shown in Figs. [5] and [6] for the standard cosmology, we present the corresponding results for the LTR cosmology in 
Figs. [7] and m 

As Figs. [5] and [7] illustrate, the dependence of fiJot on the model parameters fx, Mc, and Aq is somewhat com- 
plicated, and the results displayed therein clearly warrant detailed explanation. Perhaps the most intuitive way of 
understanding these results is to begin by examining them in certain limiting regimes. For example, consider a situ- 
ation in which Aq is relatively small. In this case, the confinement time scale to is relatively late. If to is sufficiently 
late that all modes in the tower begin oscillating immediately at to, each ^x is given in the standard cosmology by 
Eq. (|4.18|) for all A and by Eq. (|4.19|) in the LTR cosmology. In addition, let us assume that fx is large enough that 
decays can be neglected, and that Hj is at least moderately large, so that essentially all of the a\ which contribute 
meaningfully to $7*^^ survive inflation. In this special case, we can explicitly sum the contributions fix to obtain the 
result 



^It « 7T^(5gO^ T7^ 4^'Ce X , ,„ (5.6) 



feA^Y 3/21/2 J^ standard cosmology 

IbbTT \MpJ [{tG/tw)^/^ LTR cosmology . 

In deriving this result, we have used the second identity in Eq. (|2.25p . Note that Eq. (|4.14p implies that to < ^rh; 
hence il^^^ is suppressed in the LTR cosmology relative to the standard cosmology by a factor which can be quite 
significant. Indeed, this suppression factor in 51^^^ is due to the uniform suppression of each individual contribution 
il\ in this regime by the factor £ltr given for the A > l/ic case in Eq. (|4.17l) . 

Perhaps the most interesting aspect of this result is that it depends only on Aq, and is independent of both Mc 
and fx- This is quite surprising indeed, for it indicates that in this regime, no matter how many of the ax contribute 
significantly to fltot, the total contribution to the dark- matter relic abundance is the same. For example, a strongly- 
mixed scenario with y <^ 1 and a vast number of modes contribute more or less democratically to fltot will yield the 
same abundance as a weakly-mixed scenario with y 3> 1 in which a single, light axion accounts for essentially the 
entirety of the dark matter. This situation is realized in the Aq = 10 MeV and Aq — 100 MeV panels in Figs. [H] 
and [71 in which the value of i^1^^. remains essentially constant throughout the region of parameter space shown. 

At an algebraic level, the fact that Jl*^^ is a constant throughout substantial regions of {fx,Mc) space for small Aq is 
a reflection of the fact that the identity in Eq. (|2.25p holds regardless of the value of y. Of course, this identity requires 
that the sum over fix be taken over the entire KK tower, from the lowest mass eigenstate up to infinity. At a physical 
level, this is the appropriate sum to take for small Aq, because all of the modes begin oscillating at a common time to, 
and because the full structure of the tower is undisturbed by the decay of any modes which contribute meaningfully 
in the sum. Indeed, as we have seen from the panels of Figs. [5] and [3 this result is characteristic of situations in which 
A(3 is small. 

However, as we increase Aq, three effects can begin to alter this picture and thereby destroy the uniformity of fll^t'- 

• First, to becomes smaller and smaller, and consequently tx can begin to exceed to for the lower modes in the 
KK tower. In other words, these lower modes may begin to experience oscillations with staggered onset times, 
a phenomenon which begins with the lowest-lying modes in the tower and ultimately affects higher and higher 
modes as Aq increases. 

• Second, y decreases with increasing Aq (for fixed fx and Mc), and consequently more and more of the excited 
ax contribute significantly to fl*Q^■ Although the lifetimes tx of these modes also increase with increasing Aq, 
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FIG. 5: Contours of the total contribution f2*ot from a KK tower of general axions, plotted in {fx,Mc) space, assuming the 
standard cosmology. Each panel corresponds to a different choice of Ac ranging from 10 MeV to 100 TeV. In each case, we 
have taken ^ = go ~ = 1 and set Hi — 10" GeV. Note that the contour corresponding to fl^at ~ 1 has been highlighted 
with a dotted blue line in each panel for clarity. The solid oblique red line appearing in each panel indicates where y = 1, and 
proceeding from left to right, the dashed red lines correspond to y = {0.01, 0.1, 10, 100}. 
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FIG. 6: Contours of the tower fraction r;« from a KK tower of general axions, plotted in {fx , Mc) space, assuming the standard 
cosmology. As in Fig.O each panel corresponds to a different choice of Ac ranging from 10 MeV to 100 TeV; likewise, we have 
taken ^ = gc = S = 1 and set Hi = lO"'^ GeV. Once again, the solid oblique red line appearing in each panel indicates where 
y = 1, and proceeding from left to right, the dashed red lines correspond to y = {0.01, 0.1, 10, 100}. 
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FIG. 7: Same as in Fig. O but for the LTR cosmology rather than the standard cosmology. 
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FIG. 8: Same as in Fig. (6] but for the LTR cosmology rather than the standard cosmology. 
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they do so at a slower rate. As a result, a larger and larger fraction of the contributing portion of the KK tower 
is effectively truncated by decays. This effect can therefore lead to a reduction in ^2^^^, especially in the y ^ 1 
regime. 

• Third, as Ac increases, to can be pushed back into the inflationary era. The contributions from those modes 
which begin oscillating prior to or during inflation will therefore be inflated away. This too can result in a 
reduction of $1*^^ . 

Of course, which of these effects happens to be relevant in any given situation ultimately depends on the parameters 
in question, and whether we are working in the standard cosmology or an LTR cosmology. Let us therefore begin by 
examining the situation in the standard cosmology, as shown in Fig. [S] As we increase Aq from 10 MeV to 100 TeV, 
we see that a series of contours with smaller and smaller values of ^l^t emerges in the region of parameter space where 
y -^ I and gradually spreads over a substantial region of {fx, Mc) space. This is the effect of decays truncating the 
contributions from the higher modes in the tower, as discussed above. Note that for the regions of parameter space 
shown in Fig. [Sj neither of the other two effects outlined above is apparent. In particular, staggering effects only 
occur within regions of parameter space for which to < t\g , where t\^ is the time at which the lightest mode in the 
tower begins oscillating. This criterion can be rephrased as a condition on the model parameters fx, Mc, and Aq by 
substituting A^ for T in the middle line of Eq. ()4.2|) . In the y S> 1 and y ^ 1 regimes, we can approximate Xomx 
and Ao ~ Mc/2 respectively to obtain 



fx > (1.38 X 10^^ GeV) x gaC g^^'^^a) 



standard cosmology 



Mc < (8.18 X IQ-i^ GeV) x Igl^'^iAc) 



(Ao_y 



y>i 
2/<i • 



(5.7) 



Given these results, it is clear that staggering effects will not be visible in Fig. [5j for the y ^ 1 case, extremely large 
values of fx are required, regardless of the value of Aq; for the y <C 1 case, Mc < 10^'^ GeV or Aq > 100 TeV is 
required for these effects to be apparent. 

The situation is quite different in the LTR cosmology, as shown in Fig. [T] Indeed, as we increase Aq, all three of 
the above effects begin to become relevant. First, we observe the same effect of decaying a\ modes in the y ^ 1 
region that we saw in Fig. \E\ This effect is particularly evident in the Aq = 1 GeV panel of Fig. [T] However, in the 
LTR case, we also observe effects due to staggering, which begin to appear in the large- fx region. Indeed, following 
the same procedure applied above for the standard cosmology but using the top line in Eq. (|4.2p . we find that these 
effects emerge in regions of parameter space where 



LTR cosmology : < 



/_x>(1.03xlO"GcV)x5Ge 
Mc < (2.73 X 10"^^ GeV) x 



g* (Jrh) 
5*(Ag) 

3»(Ag) 



/ Trh \ I Ag \ 



gl ^{Trh) 






(5.8) 



GeV 



The first of these limiting forms accounts for the vertical contours which appear on the right side of the Ag = 100 MeV 
panel in Fig. [7] and encroach further and further to the left as Aq increases. 

Finally, as Ag grows beyond 100 GeV, we see the third effect emerging: the inflating away of heavy KK modes. 
In particular, since we have chosen Hi = 10~ GeV in this plot, the relic-abundance contributions from all modes 
with A > 3Hi/2 will be inflated away. Indeed, we see that when y ^ 1 (which implies that Aq ~ Mc/2) and 
Mc ^ 3 X 10^^ GeV, the entire tower is inflated away, yielding fij^j — 0. Indeed, since y increases with Ag for fixed 
fx and Mc, this effect spreads across a wider region as Ag increases. 

Ultimately, for large Ag, those modes which have not inflated away exhibit a completely staggered behavior in the 
LTR cosmology. This limit may be regarded as the converse of the "instantaneous turn-on" limit taken in Eq. (|5.6p 
for small Ag: indeed, we now have t\ > to for all of the modes which contribute significantly to il^^if Moreover, 
in this case the sum over fl\ can be explicitly evaluated using the first identity in Eq. (J2.25I) . allowing us to obtain 
explicit results for il^^^^ in the completely staggered limit for both the standard and LTR cosmologies: 



standard cosmology : 
LTR cosmology : 



"tot ~ 229/4^1/2 ^5gU y^j^ 



tlil.EfT^GC{y) 



2 
X , 



(5.9) 
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where C{y) = V] ;, )>}^^A\ for the standard-cosmology case. Note that the j/-dependence of C{y) is illustrated in 
Fig. 4 of Ref. [lf7 In sharp contrast with the results obtained in Eq. (|5.6p . we see that in this staggered regime, the 
expressions for 17^^^ in the standard and LTR cosmologies differ significantly. In the standard cosmology, i}^^^ depends 
non-trivially on fx, Aq, and Mj, (through its dependence on y). By contrast, in the LTR cosmology, fl*^^ depends 
on fx in this staggered-oscillation regime, but not on Aq or Mc- 

Following similar reasoning, we can also understand the behavior of the tower fraction 77,, beginning with the results 
displayed in Fig. [5] for the standard cosmology. These results nicely illustrate an important general property of 77* in 
brane/bulk models of dynamical dark matter discussed in Ref. [l| , which is that the behavior of 77* is strongly correlated 
with the value of y. In particular, when y ^ 1, the mass of the lightest mode in the KK tower is proportionally 
far lighter than those all of the excited modes, and consequently its contribution Vlx^ to the total abundance will be 
much larger than the contributions from all of those other modes combined. Indeed, this is nothing but the four- 
dimensional limit of the KK theory. This property of ry, is independent of both the specific cosmological context and 
whether or not any of the modes in the tower have staggered oscillation onset times, as can be seen from Eqs. (|4.18p 
through ()4.21|) . By contrast, when y <^ 1, each mode with a mass A <C Atrans contributes essentially equally toward 
n^Q^. when t\ = Iq for all a\, and hence 77, « 1. This behavior is manifest in the various panels of Fig. [51 in which 77, 
rapidly transitions from nearly zero to nearly unity as one crosses the y = 1 contour. 

In Fig. m we display the behavior of r]^ in the LTR cosmology. For small Ag, the situation is very similar to that 
in the standard cosmology: all modes in the tower begin oscillating at Iq, and hence 77 w for y ^ 1, while 77 « 1 
for y <^1. However, as Kq increases, the oscillation onset times for more and more of the lighter modes in the tower 
become staggered in time. In this regime, the energy densities p\ associated with the lighter modes in the tower scale 
like vacuum energy (and are hence unaffected by Hubble dilution) for a longer time before coherent oscillations set in 
and they begin to scale like massive matter, as illustrated in Fig. 1 of Ref. 1]. It follows that in this regime, the lighter 
modes account for a greater fraction of rij^j. For this reason, as discussed in Ref. [1], 77* no longer approaches unity 
for 7/ <C 1 in those regions of parameter space in which oscillation onset times are staggered, but instead asymptotes 
toi 

?7max = 1 - 4j ~ 0.189 , (5.10) 

Tr- 
ior the case of the LTR cosmology. Roughly speaking, the regions of parameter space in which this occurs are 
those in which the staggered-onset criteria in Eq. (|5.8I) are satisfied. Indeed, this effect first becomes apparent in 
the Aq = 1 GeV panel ol Fig. |S] and becomes increasingly significant as A^ increases and staggering effects become 
relevant lor smaller and smaller fx and larger and larger Mc- By Aq = 10 GeV, these staggering effects are realized 
over nearly the entirety ol (/x, Mc) space shown, leaving only a narrow strip in which all modes still begin oscillating 
at tc: and by Aq = 100 GeV, even this strip vanishes. Note also that the effect on r;* of modes being infiated away 
is apparent in the upper left of those panels in Fig. |S] for which Aq > 10 GeV. While 77 is technically undefined in 
this region ol parameter space because n^ot = 0, we have set 77, = within this region to illustrate where this effect 
is relevant. 

As discussed in the beginning ol this section, the interesting regions ol parameter space lor dynamical dark matter 
are ultimately those in which fi^^^ sa ficDM, while at the same time 77* differs significantly from zero. Given the results 
in Figs. [5] through ISl we can now determine whether this situation ever actually arises in our model. Comparing the 
results in Figs. [5]and|6l we see that this occurs in the standard cosmology lor small values ol fx, within a diagonal 
stripe ol parameter space slightly to the left ol the blue $7^^^ = 1 contour in each panel. We also see that this stripe 
moves to the right in {fx, Mc) space as Aq increases. By contrast, comparing the results in Figs. [7] and [H we see that 
the above conditions are satisfied in the LTR cosmology in the region ol parameter space where 

fx -^ 10^^ - 10^5 GeV 
preferred region (LTR) : <( • Ac > 100 GeV (5.11) 

Mc small enough that y < 1 . 



1 



Note that there are two (ultimately equivalent) ways to derive this result, corresponding to two different methods of taking the J/ — > 
limit. In Ref. [l|, we recognized that A„ ft: (n + l/2)Mc as y — > 0. Since the j; — ^ limit also implies that A <SC Atrans for all 
A, we can similarly approximate Ax ~ 1/A in this limit. We then have Qx ^ ^\ ^ l/^^i whereupon it follows that Q;^ /Htot = 
4/ J^^(n-|- l/2)~^ = 8/iT-^, or equivalently J7max = 1 — S/tt^. However, it is also possible to retain the exact form Qx ^ ^|) whereupon 
we see that fixg /f^tot = A^ where Axq is the value of Ax for the lightest eigenvalue Aq and where we have used the identity ^^ A| = 1 
to perform the sum over KK modes. Note that this result is exact and valid for all y. However, it is easy to verify that Ax„ — > 2\/2/i" 
as J/ — > 0. We thus again find that r?max = 1 — S/tt^ . 
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This result is certainly intriguing, as it suggests that the preferred scale for Aq in this model is roughly the TeV scale 
for the LTR cosmology — a scale at which there is good reason to expect new physics to appear. 

The principal message of Figs. [5] through 13 then, is that our bulk-axion model indeed satisfies the conditions on 
rijQ^ and 77, for dynamical dark matter within these regions of parameter space. In other words, within these regions, 
our axion ensemble reproduces the observed dark-matter relic abundance, and does so in a non-trivial manner, with 
a substantial number of its constituents contributing significantly to i7cDM. Of course these alone are not sufficient 
conditions for a successful model of dynamical dark matter: such a model must also not only have an appropriate 
present-day equation-of-state parameter w, , but also satisfy all additional relevant phenomenological constraints. In 
the remainder of this section, we will address the constraints on w^; the rest of the applicable constraints will be 
addressed in Sect. I VII 

One particular ramification of these constraints, however, is appropriate to mention before proceeding further. As 
discussed in Sect. IIVI certain bounds which apply generically to models with large, flat extra dimensions strongly 
prefer the LTR cosmology over the standard cosmology. For this reason, we will focus primarily on the LTR case 
from this point forward. 



C. Dark Towers: Equations of State 



Having characterized the behavior of ri^Q^. and 77* over the parameter space of our bulk-axion model, we now 
proceed to discuss the third critical quantity which characterizes the dynamical dark-matter ensemble in this model: 
the present-day effective equation-of-state parameter w^ . Since it is now clear which regions of model parameter space 
are suitable for dynamical dark matter, we will not perform a general survey of w^ over the entirety of that parameter 
space, as we did with ^*^^ and 77*, but instead focus on the preferred regions indicated in Eq. (|5.1ip . 

In order to calculate w* we need to know the values of the coefficients and exponents A, B, a, and /3 appearing in 
Eq. (15. 4p . This, in turn, requires knowledge of how our abundances and decay widths scale with A. As in the previous 
subsection, we will assume that the abundances of the a\ result from misalignment production, and likewise we will 
assume that their decay widths F^ are those appropriate for a photonic axion with c^ = 1. Because the preferred 
region of parameter space specified in Eq. (J5.11I) for our model is one which is well approximated by assuming staggered 
oscillation onset times for all relevant modes, the correct expression for ilx is the one given in Eq. (|4.2ip . Likewise, 
the decay width for a photonic axion is given by the expression in Eq. p.ip . We therefore find that the coefficients A 
and B appearing in Eq. (|5.4[) are given respectively by 
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while the power-law indices a and /3 are given by 
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Substituting these results into Eq. (15. 4p . we find that 
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Let us discuss the implications of these results. First, it was noted in Ref. [1| that the effective equation-of-state 
parameter Wcb (t) for any given dynamical dark-matter ensemble at any time t < inow will always fall within the range 
< Wcff{t) < w* as long as a + /? < — 1. This makes such ensembles less dangerous from a phenomenological point of 
view. Indeed, we see from the results above that this criterion is satisfied for both the large-A and small-A regimes in 
the bulk-axion model under consideration here. 

Second, in order to convey a sense of the characteristic size of w* in the favored region of parameter space for 
dynamical dark matter given in Eq. (|5.11|) . we note that for the choice of fx = 10^* GeV, Mc — 10^^^ GeV, and 
Ag = 1 TeV, with go = ^ ^ ^ I, we find that w* « 8.4 x IQ-^^ for A > Trmj./Mc, while w^ « 5.7 x 10"" 
for A < Trrnj^/Mc- As these numbers are both extremely close to zero, we conclude that at present time our axion 
ensemble has an effective equation of state which can be legitimately interpreted as that of dark matter. Thus our 
ensemble meets all three requirements for a self-consistent model of dynamical dark matter. 

VI. CHARACTERIZING THE ENSEMBLE: CONSTRAINTS AND PROSPECTS FOR DETECTION 

In the previous section, we demonstrated that an ensemble of mixed KK excitations of a bulk axion field can 
collectively account for the observed relic abundance of dark matter in our universe. However, as discussed in Ref. [l|, 
in order to be a viable model of dynamical dark matter, the model must also comply with a variety of additional 
laboratory, astrophysical, and cosmological constraints. Some of these constraints are intrinsic to any theory involving 
large extra dimensions, while others arise due to the physical effects of the axion field which propagates in the bulk 
of those dimensions. A number of analyses of such constraints exist in the literature [2y, |4l|, \^ for the specific case 
in which the bulk axion in question is identified with the QCD axion and the fundamental, D-dimensional quantum- 
gravity scale is taken to be roughly Md ^ ©(TeV). By contrast, in the present analysis, we are interested in a 
broader class of axions which are neither required to couple to the fields of the SM (and in particular to hadrons) 
in the same manner as a QCD axion, nor subject to the same strict relationship between the suppression scale for 
those couplings and the axion mass. Moreover, our primary motivation is not to address the hierarchy problem, but 
to address the issue of what constitutes the non-baryonic dark matter in our universe. For these reasons, we will not 
focus exclusively on scenarios in which Md is at or near the TeV scale, but also consider scenarios with much larger 
M]j. As a consequence, exclusion limits on the parameter space of the more general axion scenarios considered here 
can differ quite significantly from those presented in previous studies, and thus warrant reexamination. 

We begin our summary of the applicable constraints on our model with a brief synopsis of those limits which arise 
generically in theories with large, flat extra dimensions and which do not depend on the presence or properties of the 
bulk axion field. For the most part, these limits, an overview of which was presented in Ref. [28| . tend to derive from 
the non-observation of physical effects related to the dynamics of KK gravitons. These limits take many forms. First, 
there is the direct lower bound on Mc quoted in Eq. p.4p from experimental limits on modifications of Newton's law 
at short distances due to KK-graviton exchange [32j] . In addition, a number of constraints arise as a consequence of 
the production of these particles in the early universe [28l|36|. As discussed in Sect. lIVl these cosmological constraints 
can collectively be addressed by positing that the universe underwent a late period of cosmic infiation with a reheating 
temperature Trh ^ ©(MeV). Thus, by adopting an LTR cosmology with a reheating temperature of this order, as 
we have done, we automatically ensure that a large number of these model-independent constraints are satisfied. 

A number of additional constraints on theories of this sort can be derived from observational limits on KK-graviton 
production in astrophysical sources, such as stars [43] and supernovae |4l,|43- The most stringent of these constraints 
are currently those resulting from gravitationally trapped KK gravitons in the halos of neutron stars either decaying 
to photons or serving as a heat source for the stars themselves. In the case of n > 1 fiat extra dimensions with 
equal radii compactified on an n-torus, these limits supersede the limit on Ale given in Eq. (|3.4p . In particular, for 
n = 2, the bound is M^ > 5.8 x lO"'^ GeV, while for n = 3, one finds Mc > 3.8 x 10-i° GeV [H. However, if 
the radii of the extra dimensions differ from one another, or if the compactification manifold is not toroidal, these 
bounds can be considerably weaker. Furthermore, it is possible that the axion propagates only within some number 
Ua of the additional dimensions, Ua < n. In other words, the axion could be confined to a (4 -|- na)-dimensional 
brane within the bulk. In this case, the effective, four-dimensional scales fx and Mp are related to the fundamental, 
higher-dimensional scales fx and Md in completely different ways: 

M2 = VnMl+" 

fx = v;„/x+"" • (6-1) 

The upshot, then, is that naive limits on Md derived from KK-graviton dynamics under the assumption of toroidal 
compactification and equal radii do not necessarily translate in a straightforward manner into constraints on the mass 
scales relevant to the physics of a bulk axion. Fortunately, the bound in Eq. (13.41) is universal and is not sensitive to 
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the total number of extra dimensions, unless they are each of comparable size. We will therefore take this bound to 
be the lower limit on M^ in the ria — 1 model under consideration here. 

We now turn to address those constraints which relate to the effects of the bulk axion itself. Indeed, a number of 
considerations serve to constrain the properties of light exotic particles with suppressed couplings to SM fields. Some 
of these constraints derive from observational limits on the production of such particles in astrophysical sources such 
as stars and supernovae; others derive from limits on the decays of a cosmological population of such fields into SM 
fields; and still others owe to direct experimental bounds from microwave-cavity experiments, helioscopes, etc. A 
detailed analysis of the exclusion limits implied by these constraints on general bulk-axion scenarios will be presented 
in Ref. [35|. Here, we merely summarize the results and discuss their implications for a mixed KK tower of axions as 
a model of dynamical dark matter. 

As we shall discuss further in Ref. [35|, it is convenient to separate the applicable constraints into four rough classes, 
based on the origin of the constraint and on the dynamics being probed. The first class of constraints which apply 
to scenarios of this sort are those related to the total present-day dark- matter relic abundance ^l^t ■ Most of these 
bounds have been addressed in previous sections, but it will be useful to recapitulate them here: 

• The axion ensemble must yield an acceptable contribution to the present-day dark-matter relic density. While 
OjQ^. < OcDM is permitted, provided some additional field or fields make up the deficit, values of ^l^t i^^ excess 
of the WMAP upper bound in Eq. (|l.ll) are excluded. 

• At no time in the past may our ensemble overdose or prematurely matter-dominate the universe. 

• The present-day effective equation-of-state parameter w* for the ensemble must not deviate significantly from 
zero. 

• Misalignment production must provide the dominant contribution to 17 a for all a^, and the population of hot 
axions generated via thermal production must be negligible. We therefore require that Fpiod ^ -ff at all times 
after the end of cosmic inflation, where Fpiod is the total production rate of axions from interactions with SM 
fields in the thermal bath. 

• We have also assumed that the population of axions generated from the decays of cosmic strings associated with 
the breaking of the global U{l)x symmetry is small compared to the population generated by misalignment 
production. We therefore impose the requirement that fx ^ Hj, so that such strings are diluted away by 
inflation. 

• Our model must respect current observational limits on isocurvature fluctuations from WMAP [3| . 

The last of these constraints warrants additional discussion. Non-adiabatic fluctuations — also known as isocurva- 
ture fluctuations — refer to fluctuations not in the total energy density (which relates directly to spacetime curvature) 
but rather in how that total energy density is distributed among different contributing flelds (including the collective 
contribution from the dark sector) . Such isocurvature fluctuations are tightly constrained by a combination of CMB 
observations, baryon-acoustic-oscillation (BAO) measurements, and supernova data [3|. Such fluctuations generically 
arise whenever a cosmological population of particles is produced in a manner such that its primordial density per- 
turbations are uncorrclated with those of the inflaton fleld. Indeed, limits on isocurvature fluctuations place severe 
constraints on the relic abundance of a standard QCD axion produced via vacuum misalignment [18|, so it might 
reasonably be assumed that such limits might play a signiflcant role in constraining our model as well. 

It turns out, however, that our model satisfies the WMAP constraints on non-adiabatic fluctuations far more easily 
than do standard axion dark-matter models. A detailed discussion of these constraints and how they apply to bulk- 
axion models of dynamical dark-matter will be presented in Ref. [35{ . but the gist of the argument is as follows. 
Although our dynamical dark-matter ensemble comprises a large number of individual components a\, the fact that 
^Xat ~ ^CDM implies that the individual abundance Vt\ associated with each of these components is actually quite 
small. Furthermore, the underlying five-dimensional nature of our KK axion tower guarantees that the primordial 
density fluctuations for each a\ are all determined by the fluctuations 89 of the same initial misalignment angle 9. 
For these reasons, the expected magnitude for isocurvature fluctuations in our model turns out to be no greater than 
it is in models in which misalignment production causes a single four-dimensional field to carry the complete dark- 
matter abundance. Moreover, if one assumes a Gaussian distribution for 59, it is straightforward to demonstrate \ 
that {{69)'^) ~ H] /(2Trfx)'^- Thus, all that is required is that Hj ^ fx within our preferred regions of parameter 
space. However, as discussed in Eq. (jS.lip . the phenomenologically preferred scale for fx in our model is roughly 
0(10^'' — 10^^) GeV. For fx at or around this scale, it turns out that current constraints on isocurvature fluctuations 
can be satisfied, provided that Hj ^ 0(10^ — 10^'^) GeV. Such a scale for Hj is easy to realize in traditional cosmological 
scenarios, and is even more natural in LTR cosmologies wherein the reheating temperature is ©(MeV). Thus, in our 
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model, it is not difficult to satisfy current isocurvature bounds while simultaneously obtaining a total relic abundance 

^tot ~ ^CDM- 

The underlying reason why our model easily evades these non-adiabatic constraints is that within the preferred 
region of parameter space in Eq. (|5.1ip . the five-dimensional axion in our model is not the standard QCD axion. In 
particular, we see that the scale Aq is significantly larger than Aqcd- Our model is thus freed from the implicit para- 
metric dependence on Aqcd which afflicts more traditional models of axion dark matter, and allows the corresponding 
non-adiabatic fluctuations to have a much smaller scale. 

A second class of constraints comprises those observational limits on processes in which axions are produced via 
their interactions with the fields of the SM and then subsequently detected via those same interactions. These include: 



• 



Limits from helioscope experiments, such as CAST [46|, which search for axions produced by interactions with 
SM particles in the sun via their "conversion" to photons in the presence of a magnetic field. 



• Limits from light-shining-through- walls (LSW) experiments (see Ref. [47[ for a thorough review), including those 
by the BEV and GammaeV collaborations. 

The most stringent of these bounds is currently that from CAST; we shall therefore take the CAST bound as 
representative of this class. 

The physical processes to which this second class of limits applies are all subject to a particular effect which arises 
universally in models with both brane and bulk mass terms. This is the phenomenon of decoherence discussed in 
Refs. [l|, |26| . This decoherence phenomenon can substantially suppress the cross-sections for such processes in our 
model, and thereby significantly weaken the bounds on fx, Mc, and Aq. To summarize, the cross-section for any 
process in which axions are produced at some time to and then subsequently detected at a later time t is given by 

a{t) (X ^P(i) , (6.2) 

Jx 

where N ~ fx /Mc is the number of modes contributing in the sum and where P{t) is the detection probability at 
time t. This latter quantity is given in the relativistic limit by j26| 
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where p is the initial momentum of the axion. For any reasonable choice of model parameters, the sum in the second 
term decoheres on time scale so rapid as to be effectively instantaneous [2^1. As a result, a{t) is suppressed, relative 
to the naive expectation, by an additional factor of TV. This effect considerably weakens the constraints in this class. 
A third class of constraints can be derived from processes in which axions are produced via their interactions with 
SM fields but not subsequently detected. Instead, the presence of the axions is made manifest by their ability to carry 
away momentum and energy from a given system. These constraints include: 

• Observational limits on the energy loss in supernovae, and, in particular, on the fraction of the energy released 
by SN1987A in the form of fight exotic fields [Hi. 

• Limits related to the effects of energy dissipation by axions on stellar lifetimes. The most stringent such limits 
currently come from observations of globular-cluster stars [321 j but similar limits have also been derived from 
constraints on the lifetimes or energy- loss rates of other astrophysical bodies {e.g., the sun ^49] and white 
dwarfs fSOj'). 

• Constraints from the absence of observed signals in channels such as j + ^t and 7 -|- ^t at particle colliders. 
In general, the constraints on axion production in these channels are analogous to the well-known constraints 
on KK-graviton production [54I . 

• Limits on the branching fractions in particular exotic decay channels for certain hadrons |5l| . 

The degree to which many of these limits constrain the parameter space of bulk-axion scenarios depends quite 
crucially on how the axion in question couples to the fields of the SM. Moreover, many of the constraints in this class 
are considerably relaxed in regions of parameter space in which y ^ 1 , due to the coupling-suppression phenomenon 
discussed in Ref. [l|. This effect will be discussed in greater detail in Ref. |35{ . 

A fourth and final class of constraints is related to the interactions and decays of a cosmological population of axions. 
Depending on the cosmological epoch during which such decays occur, they can result in a number of potential signals, 
none of which have been observed to date. For example, these include: 



• 



• 
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Decays of cosmic axions which occur after the beginning of the BBN epoch (at around t ^ 1 s), but before 
last scattering (at around t ^ 10^"^ s). These could disrupt nucleosynthesis and affect the abundances of light 
elements |53j . 

Photoproduction (cither primary or secondary) from any axion decays that occur between the epoch of electron- 
positron annihilation (at around t ^ 10"^ s) and last scattering. These can lead to observable distortions of the 
CMB 0. 

Photoproduction (either primary or secondary) from any axion decays that occur after last scattering. These can 
lead to peaks and other indicative features in the diffuse X-ray and gamma-ray spectra [53 , but such features 
have not been observed by FERMI (HBl, EGRET [s^l, COMPTEL ^ [Hil, or any other X-ray or gamma-ray 
telescope [gS-il- 

• Entropy production from late axion decays. This can have observational effects on cosmological parameters, 
such as the rate of cosmic expansion. 

• Limits from microwave-cavity-detector experiments such as CARRACK [6^ and AD MX [6J| , which search for 
cosmic axions via their "conversion" to photons in the presence of strong magnetic fields. 

It should be reiterated that the vast majority of the constraints enumerated above are highly model-dependent. 
The standard energy-dissipation limit from SN1987A j48|], for example, provides one of the most stringent limits on 
the parameter space of a QCD axion. However, these limits are predicated on the assumption that the axion couples 
to hadrons with significant strength, and that processes such as NN — ?► NNa consequently dominate the axion- 
production rate. A purely photonic axion, on the other hand, lacks such couplings, and hence can only be generated 
via interactions such as the Primakoff process e~j — >■ e~a, for which the rate is much smaller. As a result, the bounds 
on fx, Mc, and Ac for such an axion are considerably weaker than those for a hadronic axion (see, for example, 
Refs. [5ll l65j for an analysis of this constraint for a four-dimensional photonic axion). A variety of other constraints, 
including bounds from monojet searches at hadron colliders and from the requirement that misalignment production 
of axions dominates over thermal production, also differ markedly depending on whether or not the axion in question 
couples to hadrons. Still other bounds, such as that from energy loss in white dwarfs [50], depend sensitively on 
whether or not a given axion couples to leptons. 

In Figs. P and [TUl we display a series of exclusion plots in {fx,Mc) space, taken from Ref. [35|, which indicate the 
regions of parameter space excluded by the considerations enumerated above. The three panels in Fig. [9] correspond 
to A(3 = {1 GeV, 1 TeV, 100 TeV} for the case of a photonic axion with Cj — I, while the three panels in Fig. [TU] 
correspond to the same choices of Ac, but for a hadronic axion with c-y — Cg — 1. In each case, we have taken 
C = .9g = ^ = 1, with Trh = 5 MeV and Hj ~ 10^'^ GeV; for the hadronic case, we have also assumed that CaTi, 
Ca-KN, etc., take the values given in Eqs. (12. 7p and (|2.8p . The shaded regions in each panel are excluded by the battery 
of constraints discussed above. The red region is excluded by CAST data, the magenta region by limits on collider 
processes in which axions appear as missing energy, the purple region by limits on modifications of Newton's law from 
Eotvos-type experiments, the orange region by limits on distinguishable features in the diffuse extragalactic X-ray 
and gamma-ray background spectra, the yellow region by observations of the lifetimes of globular-cluster stars, the 
cyan region by energy-loss limits from supernova SN1987A, the gray region by the model-consistency requirement 
that Ag < fx, and the brown region by the upper bound on the dark-matter relic abundance from WMAP. A black, 
dashed line corresponding to the condition y = n has also been included in each panel for reference. Note that each of 
the exclusion regions shown, with the exception of that from WMAP, differs from the corresponding exclusion region 
for a four-dimensional axion. The exclusion regions shown are those appropriate for the five-dimensional axion on 
which our model is based, and are derived in Ref. [Sg • 

The constraints enumerated above for which no exclusion contour has been included in these figures are generally 
subleading. For example, the applicable constraints from exotic hadron decays [5ll | are generally far weaker than the 
constraints from SN1987A, thermal production, etc., for any given choice of parameters. The constraints arising from 
observational limits on distortions of the CMB are not particularly stringent either, and turn out not to constrain any 
portion of the model parameter space shown in any of the panels appearing in Figs. [9] and 1101 This is because the 
regions of parameter space in which the T\ are sizeable are those in which fx is quite small, meaning that the Qx 
are also quite small, as is evident from Fig. [T] Constraints related to the effects of late-decaying a\ on BBN were not 
explicitly calculated in Ref. 35]. However, exclusion contours derived from BBN constraints on late-decaying particles 
are expected to be roughly similar to those derived from CMB constraints, and consequently such constraints are not 
expected to rule out any additional region of model parameter space not already excluded by other considerations. 
Limits on the effective equation-of-state parameter it;off(^) are not particularly constraining either. This should come 
as no surprise, given that we showed in Sect. |V]that the effect of decays on il^^^ was negligible within the region of 
parameter space relevant for dynamical dark matter. A number of additional constraints not listed above also serve 
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FIG. 9: Exclusion contours associated with all applicable phenomenological constraints for our bulk-axion model with Aq = 
1 GeV (left panel), Ag = 1 TeV (middle panel), and Ag = 100 TeV (right panel). In each case, we have taken ^ = go ~ i, 
with Trh ~ 5 MeV and Hi — 10^"^ GeV, and we have assumed that the axion only couples to the photon field. The shaded 
regions are respectively excluded by data from helioscope measurements with CAST (red), collider considerations (magenta), 
tests of Newton's-law modifications at Eotvos-type experiments (purple), measurements of the diffuse extragalactic X-ray and 
gamma- ray spectra (orange), observations of the lifetimes of globular-cluster stars (yellow), energy- loss limits from supernova 
SN1987A (cyan), the model-consistency requirement that Ac < fx (gray), and the upper bound on the dark-matter relic 
abundance from WMAP (brown). The black, dashed line corresponds to the condition y — n. 
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FIG. 10: Same as in Fig.[9l but for a "hadronic" axion — i.e., an axion coupled both to the photon and to the gluon field (and 
hence to pions, nucleons, etc.), but not directly to SM quarks or leptons. 



to constrain very ligfit axions and axion-like particles [66|; liowever tlie particles for whicli these constraints apply 
typically involve values of mx far smaller than those of interest here. 

It is evident from these figures that the most stringent constraints on both photonic and hadronic axions are those 
from SN1987A (cyan) and from collider limits on missing-energy processes (magenta). Nevertheless, it is also evident 
that a hadronic axion is significantly more constrained than a purely photonic axion. As discussed above, the primary 
reason for this is that the rate of axion production in a thermal setting via interactions with nuclei, pions, etc., is 
far larger than the corresponding rate of production via the electron Primakoff process and other interactions which 
involve the coupling of an axion to photons alone. 

Having assessed the phenomenological constraints on a bulk axion, we are now able to definitively address the 
question as to whether or not our model is a viable model of dynamical dark-matter. In order for this to be so, 
we require that at least some part of the preferred region in Eq. ()5.11|) be consistent with the constraints discussed 
above. Inspecting Figs. 1^ and [TUl we see that indeed our preferred region is compatible with all of these constraints in 
both the photonic and hadronic axion cases for Aq ^ 100 GeV, with Mc above the lower bound from Newton's-law 
modification, but small enough so that y ^tt. Furthermore, we also see from Figs. l9land[T0l that the phenomenological 
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constraints even permit us to reach deeply into the y <C 1 region. Note that this represents a radical departure from 
the QCD-axion results presented in Ref. [2a] — a departure which is enabled because Aq is a free parameter in our 
model. Thus, we conclude that within this region of parameter space, our bulk-axion model constitutes a viable, explicit 
model of dynamical dark matter. 

It should be stated that in addition to the limits discussed above, certain additional astrophysical bounds may also 
serve to constrain the parameter space of bulk-axion scenarios. For example, it has recently been shown [67] that 
limits on gamma-ray signals from decaying axions with masses of 0(10 — 100 MeV) produced in supernovae can yield 
an even more stringent limit than that arising from energy-dissipation considerations alone. While these bounds are 
once again mo del- dependent (and directly applicable only to cases in which the axion in question couples directly to 
hadrons with significant strength, and not to a photonic axion), they could provide an important additional constraint 
on the parameter space of dynamical dark-matter scenarios involving bulk axions. Furthermore, it is also possible 
that comparable bounds could be obtained from an analysis of photo-emission limits and cooling-rate constraints from 
neutron stars, similar to that performed for KK gravitons in Ref. |43J . 

It is important to note that while additional bounds related to axion production in supernovae may serve to further 
constrain the parameter space of bulk-axion scenarios, these constraints cannot rule out axion models of dynamical 
dark matter entirely. This is due to the fact that for any given choice of model parameters fx^ Mc, and Aq, the 
couplings of any mode for which A < nnij^/Mc to the SM fields will be suppressed by mixing effects, as discussed 
in Ref. [ll]. Indeed, because we can reach deeply into the y <C 1 region, the magnitude of this coupling suppression 
can be quite significant. For example, for y ^ 0(10"^), we find that the first twenty axion mass eigenmodes have 
coupling suppressions X^Ax ^ 10^^. If the coupling suppressions are significant for those a\ with masses in the 
"dangerous" range 10 MeV < A < 100 MeV discussed above, such a\ will be produced in supernovae at a negligible 
rate, and thus all supernova bounds on axion production can be evaded. This can be arranged by demanding that 
Trrrix/Mc > 1 GeV, so that all modes with masses A ^ 1 GeV are effectively in the small-A regime. Therefore, since 
^lat is essentially independent of Aq within our preferred region of parameter space, satisfying this condition is simply 
a matter of choosing a sufficiently large value for Aq. Indeed, setting fx — 10^^ GeV in accord with Eq. (|5.1ip . we 
find that all axion-production constraints from supernovae can be avoided for 



Ac > (56 TeV) x [j^,,r^ ) • (6.4) 

We emphasize that this rough bound is not a necessary condition for consistency with supernova data, but a sufficient 
one. Furthermore, since the neutron-star cooling and photo-emission bounds on KK gravitons rest on the assumption 
that a population of gravitationally-bound particles of this sort was generated by the supernova whose core-collapse 
produced a given neutron star, any similar bound on axions would also cease to apply in this regime. We also note that 
since the fundamental scale fx is still roughly an order of magnitude larger than the value of A^ required to satisfy 
this condition, given the input values of fx and Mc, no theoretical inconsistency results from positing a confinement 
scale of this order. 

The fact that this coupling-suppression phenomenon is capable of rendering our model consistent with supernova 
bounds despite the large multiplicity of light modes attests to the importance of this effect in brane/bulk theories. A 
more detailed overview of this phenomenon and its physical implications will be provided in Ref. [35| . 

VII. DISCUSSION AND CONCLUSIONS 

The aim of this paper has been to present an explicit realization of the dynamical dark-matter framework presented 
in Ref. [l|. To that end, we have shown that an ensemble consisting of the KK excitations of a light, axion-like field 
can indeed provide such a realization. Indeed, we have shown that despite the fact that the masses, decay widths, and 
relic abundances of all of these particles are controlled by only three dimensionful parameters, the ensemble to which 
they give rise is simultaneously able to reproduce the observed value of ficDM and satisfy all applicable constraints 
from laboratory experiments, astrophysics, and cosmology. As such, this model provides a "proof of concept" for 
dynamical dark matter as a viable alternative framework for dark-matter physics. In addition, it also provides a 
method of addressing the dark-matter question which does not require the introduction of any additional stabilizing 
symmetry. 

Many qualifications, extensions, and possible generalizations of our dynamical dark-matter framework were dis- 
cussed at the end of Ref. [l|; here, we shall restrict our attention to five points which are specific to the bulk-axion 
model presented in this paper. 



• First, in this work, we have made use of the rapid-turn-on approximation in Eq. (14.81) in calculating the relic 
abundances of the ax- As discussed in Sect. IIVI this approximation is well motivated, since the instanton- 
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generated mass term mx{T) falls rapidly with temperature when T > Ac- Furthermore, the primary results 
of this paper are essentially insensitive to this approximation. This is because the fields which contribute 
significantly to flcDM in regions of parameter space which yield a realistic dark-matter relic abundance begin 
oscillating only well after mx{T) has already settled into its constant, late-time value. However, the relic 
abundance of any field which begins oscillating before mx (T) takes this late-time value will, in general, depend 
on the details of how this mass evolves in time. The quantitative effect on the abundance of a single field has 
long been appreciated |22| . but in our model, the effects are more complicated and more subtle because we 
have a coupled system of mixed scalars with different masses and therefore different oscillation times. It would 
be interesting to examine how a more rigorous treatment of the turn-on of mx (T) would affect fitot and rj in 
situations in which these quantities are sensitive to the time-dependence of this brane-mass term. Such a study 
would have important implications for more general scenarios involving other kinds of light bulk scalars. Indeed, 
the relationship between the size of the brane-mass term for such scalars and the time at which that brane mass 
is dynamically generated may differ significantly from the relationship which holds for axions. 



• Second, as alluded to in Sect. lVIl it may be possible to further test or constrain the parameter space of bulk-axion 
models of dynamical dark matter in a number of ways. We have already mentioned one potential constraint 
which derives from limits on high-energy photons resulting from the decays of axions produced in supernovae ^63] ■ 
Other considerations also merit investigation. For example, a detailed analysis of the limits imposed by BBN 
on scenarios involving multiple decaying fields with different lifetimes and abundances could provide important 
constraints on dynamical dark-matter models in general. In addition, other considerations, such as limits on 
mass loss and decreases in the dark-matter density in the halos of dwarf galaxies 68] , could also be used to 
constrain dynamical dark-matter models. Indeed, while a number of standard constraints on individual unstable 
relic particles in the early universe have been revisited in a dynamical dark-matter context [35|, it would be 
interesting to see how other constraints would apply in this context as well. 

• Third, we note that we have not specified a particular model of inflation as part of the cosmological context for 
our model. Indeed, other than requiring a low reheating temperature Trh ^ 0{MeV), we have remained largely 
agnostic about the details of the inflationary model, the form of the inflaton potential, or even the scale Hj. 
For the most part, our model does not depend on these particulars. However, certain consistency conditions 
do place meaningful restrictions on the set of inflationary scenarios with which our model is compatible. One 
such condition can be derived from the fact that vacuum fluctuations during inflation generically give rise to a 
background value (0^) « H^ti/Air^ for any scalar (p with a mass m^ <C Hj, where ti is the duration of inflation. 
This implies that the relationship between the mass A and initial energy density px in Eq. (I4.12[) in our model is 
truly valid only for the lighter a\ in a given tower — i.e., those for which 9^A\f^ > Hftj/AiT^. By contrast, any 
heavier a\ which still satisfy A <C Hj receive the leading contributions to their background values from vacuum 
fluctuations during inflation, and thus effectively acquire an initial abundance p\ ^ X^H^tj. In typical scenarios, 
we expect Hjtj ^ N^ ^ 0(60), where N^. is the number of e-foldings of inflation. The results for iltot derived 
in Sect. FVl therefore remain consistent, provided that /|- ^ Hj. Indeed, since fx ^ 10^"* — 10^^ GeV within the 
preferred region of parameter space specified in Eq. (15. lip , we see that fx ^ Hj is certainly not inconsistent 
with our model and is in fact even expected. However, this condition on Hj has non-trivial implications for 
inflationary models. While a low scale for Hi is certainly not excluded (see, e.g., Refs. [3^, |63), extremely 
small values of Hj tend to be rather non-generic [T^ among typical classes of inflationary potentials, and thus 
require either substantial tuning or careful construction. Indeed, any consistent inflationary model of this sort 
must give rise to density fluctuations on a scale consistent with constraints from CMB data Q, such as those 
on the spectral index Us, and must also satisfy other observational constraints. The development of explicit 
inflationary scenarios of this sort is therefore an interesting topic for future investigation. 

• Fourth, we note that while we have chosen in this paper to focus on the case in which the ensemble of fields 
reproducing ricDM are the KK excitations of a bulk axion field, such a field is by no means unique in possessing 
the characteristics necessary to give rise to such an ensemble. Indeed, as discussed in Ref. [1], much of the 
analysis presented here pertains to any light bulk scalar for which a mass term is dynamically generated via 
its interactions with brane-localized fields. Furthermore, for a generic bulk scalar, the relationship between the 
time at which this mass term is dynamically generated and the magnitude of this mass term itself may differ 
from that which relates tc and mx for a bulk axion. As a result, much more freedom may exist for constructing 
viable models within the dynamical dark-matter framework. For example, light moduli could also, in principle, 
provide a viable model of dynamical dark matter. 

• Finally, we emphasize that the presence of additional axion-like fields is fairly generic, and perhaps even expected, 
in many theoretically motivated scenarios for physics beyond the Standard Model (see, e.g., Ref. ^]). Moreover, 
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it has even been argued that many of these axion-hke fields are likely to be light [66]. Thus, the discovery of a 
vast ensemble of axion-like particles could provide important insight into what physics looks like at high scales. 
Indeed, if many of these axions have relatively small masses, we find ourselves in the intriguing situation in 
which most of the matter in the universe is simultaneously both light and dark. 

Our goal in this work has been to provide an existence proof for dynamical dark matter — i.e., to provide a 
model in which lifetimes are balanced against abundances in such a way that the ensemble of dark-matter particles 
successfully reproduces ficDM while at the same time satisfying all phenomenological constraints. As we have seen in 
this paper, our bulk-axion model indeed passes this test. In one sense, our model does so in the most interesting way 
possible: with y <C 1 (signifying that our tower of axion KK modes is highly mixed) and with a tower fraction ry which 
is significantly different from zero. In another sense, however, this model is fairly conservative: those modes which 
contribute most to ^l^t turn out to be rather long-lived, and likewise our numerical result for w* within the preferred 
region of parameter space turns out to be rather close to zero. Indeed, at first glance, one might suspect that these 
latter properties are in fact generic for dynamical dark-matter models, or even that such models are therefore really 
no different from traditional dark-matter models in terms of their abundance and stability requirements. 

This is not the case, however, for the balancing of lifetimes against abundances — which is the hallmark of the 
dynamical dark-matter framework — is precisely why this framework does not require such a degree of stability, 
much less the existence of a stabilizing symmetry. While certain accidental features of our bulk-axion model result 
in a preferred region of parameter space which is somewhat conservative, we emphasize that these features are not 
generic even to theories with bulk scalars, much less realistic dynamical dark-matter models as a whole. Note, for 
example, that a particular relationship exists in bulk-axion models between the mass A of a given KK mass eigenstate 
ax, the strength of its effective coupling to SM fields, and the overall magnitude of its relic abundance i}\ through 
the dependence of these quantities on fx. Even for other bulk scalars {e.g., moduli), these relationships do not 
necessarily hold. There is therefore no reason to expect dynamical dark-matter models based around such fields to 
be as conservative as the axion model we have presented here. 

In this connection, there is an even more important point that deserves emphasis. In dynamical dark-matter 
scenarios, we have no single characteristic decay width F nor abundance Q, but rather an entire spectrum of widths 
Fa and abundances il\. This therefore begs the fundamental question: if our "proof of concept" model presented 
here is to be viewed as somewhat conservative, how far from the conservative limit can we go? 

At first glance, one might try to answer this question by attempting to determine, for each time t during the evolution 
of the universe, the maximum abundance ^muxit) that a given component in a dark-matter ensemble may have if it 
has a lifetime t ^ t. In other words, given the entirety of the cosmological constraints from BBN, CMB distortions, 
etc., there exists a function r2inax(F) which describes the maximum abundance any dark- matter constituent may have 
as a function of its decay width. It might therefore seem that knowledge of this function would uniquely determine 
the full range of possibilities inherent in our dynamical dark-matter framework. 

Such an approach to answering our fundamental question is, in a sense, already a departure from the usual manner 
of approaching dark- matter physics. However, even the notion of such a function ri,„ax(r) relies too strongly on a 
single-particle perspective. One of the critical features of our dynamical dark-matter framework is that it involves a 
vast ensemble of dark-matter components. Some of these components might decay earlier in cosmological evolution, 
while others might decay later. As a result, the maximum abundance that a given component may have if it decays on 
a characteristic time scale r will itself be directly affected not only by the abundances of all of the other components 
with earlier characteristic decay times r' < r, but even the components with r' > r. Moreover, as we have seen, 
most phenomenological constraints on dark-matter decays are sensitive not merely to what happens at a specific 
moment in time, but to the integrated effects of such decays over a broad range of time scales. In other words, our 
dynamical dark-matter framework teaches us that astrophysical and cosmological constraints do not lead to a single 
function rimax(r), but rather a more subtle set of intertwined constraints on lifetimes and abundances across our 
entire dark-matter ensemble as a whole. 

Clearly, this issue has not been studied in any detail in the literature. However, it is readily apparent that this is 
indeed the only proper way in which one should express constraints on particle decays from a generic dark sector. 
Viewed from this perspective, then, the existence of even one viable dynamical dark-matter model — no matter how 
"conservative" it might be — gives us strong motivation to re-examine cosmological and astrophysical constraints 
within this framework. Indeed, it is only in this way that we will be able to fully explore our dynamical dark-matter 
framework, and understand its full range of phenomenological possibilities. 
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Appendix A: Evolution of a Decaying Axion Field 

For completeness, in this Appendix we provide exact solutions to Eq. (|4.6I) for a real- valued function a\(t). These 
take the form 
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In these expressions, kx = \/7\ — 4A2, k was defined in Eq. (|4.7p . Af(a, 6,a;) denotes Kummer's confluent hypergeo- 
metric function 
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where (a;)„ — (x + n — l)!/(a; — 1)! is the Pochhammer function, and U{a, b, x) denotes the Tricomi confluent hyper- 
geometric function 



U{a, b, x) 



m-b) ,,,„,^, , Tib-l)^,_.^^ 



r(a-6+l) 



M{a,b,x) + 



Via) 



x'-''M{a-b+ 1,2 -b,x) 



(A4) 



It can be verified upon setting Fa = that Eq. (|A1[) reduces to the exact form obtained for a tower of stable KK 
axions in Ref. [26]. 

The values oi c]^ and c-^ in Eq. (jAip are determined by the initial conditions chosen for ax{t) and dx(t) at i = tg, 
where tg is some initial time. Expressed in terms of these initial values, this equation takes the general form 
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, (A5) 
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where the time derivatives of Af„(A,i) and U^{X,t) have the explicit forms 
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Once again, if we set Fa = in this expression (which also implies that fcA = 2iA), we recover the result 
axit) 
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which agrees with the result obtained in Ref. \2&\. 

In the rapid-turn-on approximation, in which mx{t) takes the Heaviside form specified in Eq. (j4.8p . the initial 
conditions for a\ and a\ at to = t\ take the form given in Eq. (j4.9|) . Upon substituting these initial conditions into 
Eq. (jA5[) . we find that during the cosmological epoch in which coherent oscillations of a given a\ begin, we have 

The value of any ax during subsequent epochs can then be obtained iteratively from this relation. 
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